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SECTION  1 


INTRODUCTION 

The  design  of  components  such  as  nozzles  and  reentry  nosetips  which  are 
subjected  to  a  high  pressure  and  severe  thermochemical  environment  requires  a 
knowledge  of  the  environment  and  the  properties  of  the  materials  to  be  used. 

In  addition,  it  requires  a  set  of  tools  which  can  use  this  information  to 
predict  the  thermochemical  and  thermostructural  performance  of  the  particular 
component.  From  a  technical  standpoint,  the  sophistication  which  is  required 
in  the  design  and,  in  particular,  in  the  capabilities  of  the  thermochumical  and 
thermostructural  design  tools  depends  on  the  performance  requirements  of  the 
component  in  relation  to  the  capabilities  of  the  available  materials.  If  a 
simplified  analysis  shows  large  margins  of  safety,  then  a  more  sophisticated 
.approach  is  not  warranted.  On  the  other  hand,  if  small  or  negative  margins  of 
safety  are  predicted  by  the  simplified  approach,  a  more  detailed  investigation 
is  warranted.  Since  sophisticated  analyses  are  generally  time  consuming  and 
expensive,  the  cost  of  the  analysis  must  be  weighed  against  the  hardware  and 
test  costs.  If  these  costs  are  cheap  relative  to  the  analysis  costs  and  if  the 
schedule  allows  a  cut  and  try  approach,  then  performing  a  sophisticated  analysis 
often  times  can  not  be  justified.  However,  if  the  reverse  is  true,  then  the 
more  sophisticated  approach  to  the  design  problem  requires  a  detailed  analysis. 

One  of  the  more  difficult  design  problems  which  requires  this  more  sophisti¬ 
cated  approach  is  the  use  of  pyrolytic  graphite  (PG)  coatings  in  the  throat 
region  of  either  large  diameter  solid  rocket  nozzles  which  are  intended  to 
operate  at  moderate  cheimber  pressures  or  small  dieuneter  nozzles  which  are  in¬ 
tended  to  operate  at  high  chamber  pressures.  Existing  analytical  tools  in  the 
form  of  computer  programs  have  been  validated  for  predicting  the  thermochemical 
performance  of  these  inserts.  These  predictions  include  the  definition  of  the 
surface  recession  rates  and  the  in-depth  temperature  distributions.  However, 
the  prediction  of  the  thermostructural  response  of  the  material  is  not  as  far 
advanced.  Sophisticated  analysis  procedures  are  available  in  the  form  of  computer 
programs,  but  for  many  materials  or  combinations  of  materials,  these  analysis 

V 

procedures  have  not  been  validated.  The  problem  of  validating  the  analysis 
procedures  is  generally  two-fold.  The  first  part  is  the  lack  of  meaningful 
experimental  data  for  comparison  and  the  second  part  is  the  lack  of  accurate 
high  temperature  mechanical  properties. 


Even  with  the  above  mentioned  difficulties/  a  sophisticated  analysis 
using  the  best  available  properties  can  provide  insight  into  the  thermostructural 
problem  by  identifying  what  effects  are  important.  In  this  regard,  a  detailed 
thermostructural  analysis  was  performed  for  a  seven  inch  diameter  pyrolytic 
graphite  coating  on  an  AGSR  substrate.  The  items  investigated  Included  the 
effect  of  grooves  in  the  substrate  on  the  cooldown  manufacturing  stresses, 
the  effect  of  performing  an  elastic-plastic  analysis  versus  a  linear  elastic 
analysis,  and  the  effect  of  backside  constraint  on  the  stress  distribution  in 
the  insert  during  the  motor  firing.  Another  significant  result  was  identifying 
the  validity  of  presently  available  property  data. 

The  significant  results  and  conclusion  obtained  from  tne  thermochemical 
and  thermostructural  analyses  are  summarized  in  Section  2.  The  material  property 
data  used  in  the  analysis  are  presented  and  discussed  in  Section  3.  The  thermo¬ 
chemical  results  presented  in  Section  4  are  followed  by  the  thermostructural 
results  in  Section  5.  The  thermostructural  results  are  evaluated  in  Section  6. 
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SECTION  2 


SUMMARY  OF  RESULTS  AND  CONCLUSIONS 

A  thermostructural  analysis  was  performed  for  a  7  inch  diameter  PG  coated 
nozzle  throat  insert.  The  plan  for  this  analysis  is  summarized  in  Figure  2-1. 

The  conclusions  which  resulted  from  the  analyses  outlined  in  this  figure  are 

•  The  effect  of  heat  conduction  through  the  upstream  and  downstream 
bulk  graphites  has  a  significant  influence  on  the  AGSR  temperature 
and  thus  on  the  strain  distribution  in  the  insert  but  a  negligible 
effect  on  the  PG  coating  surface  recession  rate 

•  The  lip  on  the  upstream  bulk  graphite  is  marginal  for  the  protection 
of  the  leading  edge  of  the  PG  coating 

•  The  inclusion  of  the  PG  surface  recession  rate  in  the  structural  analy¬ 
sis  would  be  expected  to  have  a  significant  influence  on  the  late 

time  predicted  insert  strain  distributions 

•  The  inclusion  of  grooves  in  the  substrate  has  a  beneficial  effect  on 
the  PG  interlaminar  shear  stress  but  a  detrimental  effect  on  the 
stresses  in  the  AGSR 

•  The  elastic-plastic  analysis  showed  significant  regions  of  plastic 
flow  in  the  AGSR  and  predicted  a  larger  region  of  failure  in  the  AGSR 
than  the  liner  elastic  analysis 

•  Predicted  modes  of  failure 

-  Tensile  failure  of  AGSR  during  cooldown 

-  Compressive  failure  of  AGSR  at  a  motor  burn  time  of  1.5  seconds 

-  Compressive  failure  of  AGSR  and  PG  coating  at  a  motor  burn  time  of 
60  seconds  (5  mil  displacement) 

•  For  constant  backside  constraint,  the  motor  firing  stresses  are  a 
maximum  just  prior  to  motor  burnout 

•  Superposition  of  elastic  solutions  is  not  valid  for  materials  which 
have  temperature  dependent  properties 

•  It  is  probable  that  the  data  banks  used  in  the  analysis  for  both  the 
AGSR  and  the  PG  coating  are  inconsistent  (i.e.,  not  based  on  a  con¬ 
sistent  set  of  material  samples) 
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In  order  to  improve  the  thermostructural  analysis  procedure,  effort  must 
be  initiated  to 

•  Characterize  the  materials  being  used  and  understand  the  effect  of 
varying  the  primary  fabrication  process  variables 

s  Perform  an  Integrated  analysis  and  test  program  to  obtain  experimental 
data  which  can  be  used  to  validate  the  design  approach. 
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SECTION  3 


'  NOZZLE  CONFIGURATION 

The  7.0  inch  diameter  nozzle  is  a  submerged  nozzle  composed  of  ablative 
materials  which  were  selected  based  on  both  their  thermal  insulating  and  ablative 
characteristics.  The  nozzle  geometry  and  materials  are  described  briefly  in 
Section  3.1,  and  the  thermal  and  mechanical  properties  required  to  perform  a 
thermostructural  analysis  are  presented  in  Section  3.2. 

3.1  NOZZLE  GEOMETRY 

The  geometry  and  materials  for  the  7-inch  diameter  submerged  nozzle  are 
shown  in  Figure  3-1.  In  this  nozzle  design,  rubber  insulation  and  MX  2600  silica 
phenolic  are  used  on  the  backside  of  the  nozzle  where  the  gas  velocities  are 
relatively  low.  The  nose  cap  is  fabricated  from  MXCn  280  vhich  is  an  elastomer 
modified  carbon  phenolic.  Bulk  graphite  grades  Speer  8882  and  9139  protect  the 
upstream  and  downstream  ends  of  the  nozzle  throat  insert  which  is  composed  of 
a  pyrolytic  graphite  coating  (PG)  on  an  AGSR  substrate.  The  throat  insert  and 
the  bulk  graphites  are  insulated  from  the  steel  structural  support  by  MX  2600 
silica  cloth  phenolic.  A  throat  exit  ring  of  MX  4926  carbon  cloth  phenolic 
at  a  layup  angle  of  90®  is  located  downstream  of  the  Speer  9139  graphite. 

The  details  of  the  nozzle  throat  assembly  are  shown  in  Figure  3-2  and  the 
details  of  the  PG  coated  nozzle  throat  insert  are  shown  in  Figure  3-3.  The  more 
important  assembly  details  shown  in  Figure  3-2  are  1)  the  AGSR  is  bonded  to 
both  of  the  bulk  graphites  using  Armstrong  A-2  epoxy  resin  and  using  a  bondline 
thickness  from  5  to  20  mils,  2)  the  AGSR  substrate  is  overwrapped  with  3  layers 
of  2  mil  thick  teflon  tape  and  3)  an  8  to  10  mil  thick  bondline  of  A-2  epoxy 
resin  is  used  between  the  teflon  overwrap  and  the  MX  2600  silica  phenolic  insula¬ 
tion.  In  Figure  3-3,  the  details  of  the  PG  coated  insert  prior  to  final  machining 
are  shown.  Grooves  are  located  in  the  upstream  and  downstream  ends  of  the  AGSR 
substrate.  The  structural  analyses  presented  in  Section  5  investigate  the  effects 
of  these  grooves  on  the  resulting  stress  distribution  in  the  throat  insert.  As 
shown  in  Figure  3-3,  the  PG  coating  has  a  thickness  of  60  to  65  mils  upstream 
of  the  nozzle  throat  but  this  thickness  decreases  linearly  to  45  to  50  mils 
between  the  throat  and  insert  exit.  Prior  to  installing  the  throat  insert  into 
the  nozzle  assembly,  the  final  step  in  its  fabrication  is  the  cleanup  machining 
on  both  ends  of  the  insert. 


Two  items  to  note  in  Figure  3-1  relative  to  the  flow  field  and  boundary 
layer  analyses  are;  1)  the  boundary  layer  initiation  is  assumed  to  occur  at  the 
interface  between  the  silica  phenolic  and  the  MXCE  280  carbon  phenolic  and 
2)  the  step  downstreeun  of  the  bulk  graphite  throat  Insert  entrance  ring  is  modi¬ 
fied  by  the  dashed  line  shown  in  this  figure.  The  effect  of  the  boundary  layer 
initiation  location  on  the  thermal  and  structural  results  is  relatively  minor. 

The  second  assumption  is  justifiable  based  on  the  fact  that  the  objective  of 
the  analysis  is  to  obtain  an  overall  assessment  of  the  thermal  and  structural 
capabilities  of  the  nozzle  throat  insert  and  is  not  an  evaluation  of  a  particular 
local  perturbation.  This  local  perturbation  would  be  of  primary  significance 
in  predicting  the  local  surface  recession  rate  and  not  in  predicting  the  in- 
depth  transient  temperatures. 

3 . 2  MATERIAL  PROPERTIES 

The  material  properties  required  to  perform  both  the  chemical  and  struc¬ 
tural  analyses  are  presented  in  this  section.  A  thermal  analysis  was  performed 
for  both  the  nozzle  throat  assen^ly  (Figure  3-2)  and  the  nozzle  throat  insert 
(Figure  3-3)  .  Therefore,  thermal  properties  are  required  for  the  graphititic 
and  ablative  materials  shown  in  Figure  3-2.  These  properties  are  presented  in 
Section  3.2.1.  The  structural  analysis  was  performed  only  for  the  nozzle  throat 
insert,  and  thus,  mechanical  and  thermal  expansion  properties  are  required  only 
for  the  PG  coating  and  the  AGSR  substrate.  These  properties  are  presented  in 
Section  3.2.2.  One  of  the  parameters  studied  in  the  structural  analysis  was 
the  constraint  on  the  backside  of  the  AGSR.  In  this  study,  a  thin  cylindrical 
shell  of  multi-layered  grafoil  was  considered  in  order  to  obtain  a  flexible 
backside  constraint.  The  properties  used  for  this  material  are  also  presented 
in  Section  3.2.2. 

3.2.1  Material  Thermal  Properties 

The  materials  considered  in  the  two  dimensional  heat  conduction  thermal 
analysis  are  the  Speer  8882,  Speer  9139,  and  AGSR  bulk  graphites,  the  PG 
coating,  and  the  MX  2600  silica  and  MX  4926  carbon  cloth  phenolics.  The  pro¬ 
perties  required  in  the  thermal  analysis  arc 

•  thermal  conductivity  in  both  directions  as  a  function  of  temperature 

•  specific  heat  as  a  function  of  temperature 

•  density 

The  thermal  conductivity  of  the  PG  coating  as  a  function  of  temperature 
is  shown  in  Figure  3-4  for  both  the  "a"  and  "c"  directions.  Two  different  curves 
are  shown  in  this  figure  for  the  "c"  direction  and  the  reference  for  each  set 
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of  data  is  noted.  The  difference  between  the  two  dets  of  "c"  direction  thermal 
conductivities  is  significant,  and  the  effect  of  this  difference  on  the  predicted 
in-depth  transient  temperatures  is  presented  in  Section  4.  In  the  text  that 
follows,  the  curve  designated  by  the  solid  line  in  Figure  3-4  is  referred  to  as 
the  revised  "c”  direction  thermal  conductivity,  and  the  dashed  line  is  referred 
to  as  the  original  "c”  direction  thermal  conductivity.  The  thermal  conductivity 
values  of  the  bulk  graphites  and  charring  ablating  phenollcs  are  presented  in 
Figures  3-5  and  3-6,  respectively.  The  values  for  the  Speer  graphites  were  ob¬ 
tained  from  Reference  3  and  the  nominal  values  for  the  AGSR  were  obtained  from 
Reference  4.  For  the  majority  of  the  temperature  range,  the  with  and  across 
grain  values  of  the  AGSR  thermal  conductivity  were  obtained  by  applying  a  i  20 
percent  variation  to  the  nominal  values.  For  the  cloth  phenollcs,  the  curves 
in  Figure obtained  |rom  Reference  5  and  represent  a  combination  of  the 
virgin  yilhd  char  thermal  cond^ictlvities. 

The  values  used  for  the  specific  heat  of  these  materials  are  presented 
in  Table  3-1.  The  specific  heat  values  for  the  PG  coating  were  obtained  from 
Reference  2;  the  specific  heats  of  the  AGSR  and  Speer  graphites  were  assumed  to 
be  the  same  as  ATJ;  and  the  specific  heats  for  cloth  phenollcs  were  obtained 
from  Reference  5.  For  these  latter  materials,  the  values  for  specific  heat  in¬ 
clude  the  effect  of  the  energy  eJssorbed  by  the  in-depth  material  degradation. 

The  reason  for  including  this  effect  in  the  thermal  properties  is  that  the  two 
dimensional  heat  conduction  code  used  in  the  thermal  analysis  does  not  have  the 
option  of  in-depth  material  degradation.  As  a  result,  effective  thermal  pro¬ 
perties  which  include  the  effect  of  material  degradation  must  be  determined  from 
the  virgin  and  char  thermal  properties  of  the  ablative  materials.  The  den¬ 
sities  for  all  the  materials  considered  in  the  thermal  analysis  are  presented 
in  TeUsle  3-2. 

3.2.2  Structural  Material  Properties 

The  structural  analyses  to  be  described  in  Section  5  for  the  nozzle  in¬ 
sert  composed  of  the  PG  coating  and  the  AGSR  substrate  were  performed  considering 
both  elastic  and  plastic  effects.  The  mechanical  properties  used  in  the  elastic 
structural  analysis  are  presented  in  Section  3. 2. 2.1,  and  the  additional  mechani¬ 
cal  properties  required  for  the  plastic  structural  analysis  are  presented  in 
Section  3. 2. 2. 2.  Some  guidelines  concerning  the  accuracy  of  the  material  pro¬ 
perties  used  in  the  structural  analysis  are  presented  in  Section  3.2.2. 3,  and 
information  concerning  the  strength  of  the  PG  and  AGSR  is  presented  in  Section 
3. 2. 2. 4. 
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3. 2. 2.1  Elastic  Structural  Analysis  Properties 

For  an  axisymmetric  elastic  structural  analysis,  the  constitutive  equa¬ 
tions  in  cylindrical  coordinates  for  orthotropic  materials  are 


where  E^,  Eg  -  Young's  molulus 

-  shear  modulus 
rz 

''zr'  ''rz'  ''er'  '’rO'  '’ze  “  Poisson’s  coefficients 

a..*  OL  ,  a_  -  coefficients  of  linear  thermal  expansion 

If  in  addition  to  being  orthotropic,  the  materials  have  a  plane  of  isotropy 
then  they  are  referred  to  as  transversely  isotropic  and  the  following  relations 
exist  between  the  constants  in  the  stress-strain  constitutive  relations  described 
by  Eq  (3-1)  . 


Sj,  for  r,  9  plane  of  isotropy  '’zr 


\)  V 

zr  rz 

z  r 


E  for  z,  0  plane  of  isotropy 
z 


(3-2) 


or  E  *  E  for  r,  z  plane  of  isotropy 
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Equations  (3-1)  and  (3-2)  have  five  independent  stress-strain  constants  and  3 
independent  coefficients  of  thermal  expansion,  and  these  constants  are  required 
as  input  to  the  structural  analysis.  For  a  cylindrical  tube  of  AGSR  coated 
with  PG,  (refer  Figure  3-7a) ,  the  plane  of  isotropy  for  the  PG  is  the  2,6  plane 
and  the  plane  of  'sotropy  for  the  AGSR  substrate  is  the  r,6  plane.  For  the  PG 
coating,  the  z  and  6  directions  are  the  with  grain  directions  (often  referred 
to  as  the  "a"  direction)  and  the  r  direction  is  the  across  grain  direction 
(often  referred  to  as  the  "c"  direction).  For  the  AGSR,  the  r  and  6  directions 
are  the  across  grain  directions  and  the  z  direction  is  the  with  grain  direction. 
For  the  billet  size  required  for  the  7  inch  diameter  insert,  this  material  is 
generally  considered  isotropic.  For  the  actual  nozzle  insert,  the  coordinate 
axes  used  in  the  analysis  is  shown  in  Figure  3-7b.  The  plane  of  isotropy  for 
the  AGSR  in  this  configuration  is  the  same  as  described  above.  However,  the 
plane  of  isotropy  for  the  PG  coating  is  the  2-0  plane  which  is  a  curved  sur¬ 
face.  The  1-direction  in  this  coordinate  system  corresponds  to  the  "c"  direc¬ 
tion  of  the  PG  coating  and  the  2  and  0  directions  correspond  to  the  "a"  direc¬ 
tion  of  the  coating.  In  using  the  DOASIS  computer  program,  the  "a”  and  "c" 
direction  material  properties  are  input,  and  the  program  rotates  them  to  the 
r  and  z  directions 

The  AGSR  and  PG  coating  thermal  expansion  ilata  used  in  the  analysis  arc 
presented  in  Figure  3-8  as  a  function  of  temperature  for  both  the  with  and  across 
grain  directions.  The  PG  coating  data  were  obtained  from  Reference  1.  How¬ 
ever,  Reference  1  obtained  the  data  from  Referencr-  6  and  it  appears  as  though 
Reference  6  obtained  it  from  Reference  7.  The  AGSP  data  were  also  obtained 
from  Reference  1.  These  data  for  the  AGSR  were  determined  based  on  known  room 
temperature  values  and  on  assuming  that  the  temperature  dependency  was  the  same 
as  for  similar  grades  of  graphites. 

Young's  moduli  as  a  function  of  temperature  for  the  PG  are  shown 
in  Figure  3-9  for  both  the  with  and  across  grain  directions.  In  addition,  this 
figure  shows  the  dependence  of  the  "a"  direction  Young's  modulus  on  the  direc¬ 
tion  of  the  applied  stress.  The  Young's  modulus  values  in  the  "c"  direction 
and  the  tension  Young's  modulus  values  in  the  "a*  direction  were  the  same  as 
used  in  Reference  1  and  were  obtained  from  Reference  8.  The  temperature  depen¬ 
dency  of  the  "c"  direction  modulus  was  based  on  extrapolating  the  measured  room 
temperature  value  by  the  ratio  (E^/E^) .  The  "a"  direction  compression  modulus 
values  tfere  obtained  frexn  Reference  9  with  the  original  source  of  the  data 
being  Reference  10. 


Three  sets  of  data  for  the  Young’s  moduli  of  AGSR  graphite  are  shown 
in  Figure  3-10.  The  set  ihown  by  the  solid  lines  were  obtained  from  Reference 
1  and  reflect  considerable  material  anisotropy.  The  temperature  variation 
of  this  data  was  obtained  by  extrapolating  room  temperature  values  based  on 
measurements  on  similar  grades  of  graphite.  The  data  points  shown  in  Figure 
3-10  reflect  uniaxial  tensile  measurements  on  cross  grain  AGSR  samples  taken 
from  24  inch  diameter  billets.  Room  temperature  data  presented  in  Reference  12 
indicate  that  the  AGSR  is  essentially  isotropic  for  billets  in  the  14  to  35 
inch  diameter  range.  The  data  shown  in  Figure  3-10  do  not  reflect  any  varia¬ 
tion  of  the  Young's  modulus  with  temperature  for  temperatures  less  than  4000°F. 
The  mechanical  properties  presented  in  Reference  11  shows  this  same  trend  with 
temperature.  The  values  used  in  the  structural  analysis  were  obtained  from  the 
third  curve  presented  in  Figure  3-10  and  were  assumed  to  be  independent  of  grain 
direction.  This  curve  was  obtained  from  the  values  represented  by  the  dashed 
curve  through  the  data  by  multiplying  these  values  by  a  factor  of  0.8.  The 
purpose  of  this  factor  was  to  account  for  the  specimens  being  in  an  unannealed 
state  whereas  the  substrate  used  in  the  nozzle  was  annealed. 

The  room  temperature  Poisson's  ratios  used  in  the  analysis  for  the  PG 
and  AGSR  are  listed  in  Table  3-3,  and  because  no  data  were  available  showing 
the  functional  relationship  with  temperarure,  they  were  assumed  to  be  independent 
of  temperature.  The  values  for  the  PG  were  obtained  from  Reference  10  and  those 
for  the  AGSR  were  obtained  from  Reference  12. 

The  remaining  constants  to  be  defined  are  the  shear  modulus  (or  modulus 
of  rigidity)  for  both  the  PG  coating  and  the  AGSR.  Since  no  data  were  available, 
the  s.'tear  modulus  for  the  AGSR  was  evaluated  from  the  relationship 


G 


‘  rn&rr 


(3-2) 


For  the  PG  coating,  a  value  of  0.2  x  10*  psi  was  obtained  from  Reference  13. 
This  value  is  in  general  agreement  with  the  one  (0.1  x  10*  psi)  reported  in 
Reference  9  which  was  obtained  by  assuming  that  the  difference  between  F  in 
pure  tension  and  E  in  beam  deflection  is  due  to  a  shear  component  in  the  beam 
deflection  method. 


In  one  of  the  structural  analyses,  a  thin  piece  of  multi-layered  grafoil 
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was  placed  on  the  outside  surface  of  the  AGSR  suooLrate.  The  properties  as¬ 
sumed  for  this  material  are  presented  in  Table  3-4. 

3. 2. 2. 2  Plastic  Structural  Analysis  Properties 

The  plasticity  model  used  in  the  DOASIS  structural  analysis  code  is  based 
on  the  total  strain  or  deformation  theory  and  more  information  concerning  the 
theory  and  its  implementation  in  the  DOASIS  code  can  be  found  in  References 
14  through  16.  For  this  theory,  the  total  strains  are  defined  as 

E  P 

e  =  (3-3) 


where 

E 

E  -  total  elastic  strain 

p 

e  -  total  plastic  strain 

The  total  elastic  strains  are  defined  by  the  set  of  constitutive  equations  given 
in  Eq  (3-1).  The  equations  defining  the  total  plastic  strains  can  be  re¬ 
presented  in  matrix  form  by  the  relationship 

-P 

^  [a]  o  (3-4) 

0 

where 

-P 

c  -  effective  plastic  strain 

0  -  effective  stress 

[a]  -  matrix  of  anisotropic  parameters 

0  -  stress  vector 

For  a  stress-strain  curve  which  can  be  approximated  by  a  bilinear  fit  such  as 
shown  in  Figure  3-11,  the  mechanical  properties  in  addition  to  those  presented 
in  Section  3.2.2  which  are  required  to  solve  Eq  3-4  for  the  total  plastic  strains 
are  the  tangent  moduli  and  yield  stresses  of  the  material.  As  shown  in  Figure 
3-11,  the  tangent  modulus  is  determined  by  the  strain  hardening  portion  of  the 
stress  strain  curve  and  the  yield  stress  is  determined  by  the  intersection  of 
the  linear  segments  of  the  bilinear  approximation  of  the  stress  strain  curve. 

A  least  squares  curve  fit  procedure  is  used  to  obtain  this  bilinear  approximation. 

Since  the  plastic  structural  analysis  is  performed  for  only  the 
cool  down  of  the  insert  following  the  PG  coating  operation,  the  tangent  moduli 
and  the  yield  stresses  are  required  only  at  room  temperature.  However,  where 


experimental  stress-strain  data  are  available  at  other  temperatures,  these 
data  are  used  to  define  the  variation  of  the  tangent  moduli  with  tempera¬ 
ture. 

For  the  AGSR  graphite,  the  uniaxial  stress-strain  curves  presented  in 
Reference  11  were  used  as  the  basis  for  defining  the  bilinear  curve.  These 
data  and  the  corresponding  bilinear  fits  are  presented  in  Figures  3-12,  3-13, 
and  3-14  for  specimen  temperatures  corresponding  to  70,2000,  and  3425®F.  The 
Young's  and  tangent  moduli  obtained  from  these  curve  fits  are  summarized  in 
Figure  3-15.  No  experimental  data  were  found  for  defining  the  shear  tangent 
modulus  so  this  value  was  estimated  based  on  the  ratio  of  the  Young's  and 
tangent  moduli. 

For  the  PG,  the  uniaxial  stress  strain  curves  presented  in  Reference  6 
were  used  as  the  basis  lor  defining  the  required  mechanical  properties.  These 
curves  and  their  corresponding  bilinear  curve  fits  for  the  "a"  direction  and 
for  various  specimen  temperatures  are  presented  in  Figure  3-16  for  tensile 
stresses  and  in  Figure  3-17  for  compressive  stresses.  The  Young's  and  tangent 
moduli  obtained  frort  these  bilinear  curve  fits  are  presented  in  Figure  3-18. 

For  the  room  temperature  values  of  the  "c”  direction  and  shear  tangent  moduli, 
no  experimental  dai-.a  were  available  so  these  values  of  the  tangent  moduli  were 
obtained  using  the  ratio  of  the  "a"  direction  tensile  Young's  and  tangent 
moduli. 

The  room  temperature  yield  stresses  used  in  the  analysis  for  the  AGSR 
and  PG  coating  are  summarized  in  Table  3-5.  As  mentioned  previously,  these 
yield  stresses  correspond  to  the  intersection  of  the  linear  segments  in  the  bi¬ 
linear  stress-strain  curve  and  thus  can  be  influenced  by  the  method  used  to  per¬ 
form  the  curve  fitting. 

3. 2. 2. 3  Evaluation  of  Material  Properties 

It  is  obvious  from  Sections  3.2.1,  3. 2. 2.1,  and  3. 2. 2. 2  that  the  material 
properties  for  both  the  AGSR  graphite  and  the  PG  were  not  obtained  from  a 
single  source  but  rather  were  gathered  from  a  number  of  sources.  Because  of 
this,  it  is  pertinent  to  review  some  of  the  data  presented  in  these  sections 
and  to  make  an  evaluation  of  the  variations  that  might  be  expected  in  the 
material  properties.  These  variations  are  presented  first  for  PG  and  then  for 
the  AGSR. 

A  review  of  the  fabrication  process  of  the  PG  coating  revealed  that 
several  variables  (deposition  temperature  and  pressure,  gas  flow  rate,  type 
of  substrate  surface  finish,  coating  thickness,  deposition  of  soot  particles, 
etc)  are  significant  in  determining  the  material  properties  of  the  end  product. 
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However,  the  variation  in  material  properties  presented  here  is  restricted 
primarly  to  the  differences  represented  by  substrate  and  continuously  nucleated 
types  of  PG.  Comparison  of  existing  data  for  these  two  types  of  PG  are  shown 
in  Figure  3-19  for  the  "a"  direction  thermal  expansion  curves  and  in  Figure 
3-20  for  the  "c"  direction.  One  of  the  curves  in  each  of  these  figures  was  taken 
from  the  data  presented  in  Section  3. 2. 2.1,  the  substrate  nucleated  and  the  con¬ 
tinuously  nucleated  curves  were  taken  from  Reference  17,  and  the  data  showing  the 
variation  from  the  top  to  the  bottom  layer  in  the  PG  coating  were  taken  from 
Reference  9.  These  comparisons  show  that 

•  The  type  of  coating  nucleation  can  have  a  20  percent  effect  on  the 
thermal  expansion  of  the  PG 

•  The  location  of  the  PG  (i.e.,  whether  it  is  close  to  t^  mandrel  or 

substrate  or  close  to  the  surface)  can  have  a  15  percent  effect  on 
the  thermal  expansion  of  the  PG  ^ 

•  The  reported  thermal  expansion  data  differ  by  as  much  as  a  factor  of 
two  for  "a"  direction  but  by  only  15  percent  for  the  "c"  direction 

Since  the  mismatch  of  the  "a”  direction  thermal  expansion  between  the  coating 
and  the  substrate  is  the  dominate  cause  of  the  hoop  stress  in  the  coating,  the 
possible  variation  of  this  property  between  different  types  of  PG  is  particularly 
significant. 

Figure  3-21  shows  the  variation  of  the  Young's  modulus  for  different  types 
of  coatings.  The  variation  in  this  property  for  just  the  tension  values  is  on 
the  order  of  75  percent,  and  this  variation  occurs  primarily  at  low  and  high 
material  temperatures.  In  addition  to  the  variation  in  material  properties 
mentioned  above.  Reference  18  indicates  that  fabrication  process  variables  can 
cause  a  factor  of  two  variation  in  the  "a"  and  "c"  direction  thermal  conductivities 
and  Reference  17  indicates  fabrication  process  variables  can  cause  a  factor  of 
two  variation  in  the  torsional  strength  of  the  PG  coating. 

After  the  structural  analyses  were  well  underway,  the  data  presented  in 
Reference  19  were  obtained  for  AGSR.  These  data  included  stress-strain  curves 
from  shear  and  tensile  specimens  in  the  with  and  across  grain  direction.  The 
with  and  across  grain  curves  are  compared  in  Figures  3-22  with  the  curve  pre¬ 
sented  previously  in  Figure  3-12.  The  primary  difference  between  the  data  pre¬ 
sented  in  Reference  9  and  the  data  used  in  the  analysis  is  that  the  Young's 
modulus  used  in  the  analysis  was  assumed  the  same  in  the  with  and  across  grain 
directions  whereas  the  data  from  Reference  16  shows  a  significant  material  aniso¬ 
tropy  in  the  different  grain  directions. 
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3. 2. 2.4  Mechanical  Strength  Properties 

The  mechanical  property  data  presented  In  Sections  3. 2. 2. 2  and  3. 2. 2. 3 
are  sufficient  for  predicting  the  stress  or  strain  distribution  in  the  nozzle 
insert.  However,  except  for  the  stress-strain  curves  which  show  specimen  frac¬ 
ture,  these  data  are  not  sufficient  for  predicting  material  failure.  In  reveiw- 
ing  the  literature,  the  majority  of  the  strength  data  for  PG  is  presented  in 
the  form  of  an  ultimate  stress  limitation  based  on  data  from  uniaxial  stressed 
specimens.  Since  the  stress-strain  curves  for  most  graphites  are  nonlinear,  the 
ultimate  strain  provides  a  more  realistic  basis  for  evaluating  the  structural 
capability  of  the  insert.  This  is  particularly  true  when  evaluating  the  results 
of  an  elastic-plastic  analysis.  The  following  paragraphs  present  the  ultimate 
stress  and  strain  data  which  were  found  in  the  literature  and  which  were  used 
to  evaluate  the  structural  results  presented  in  Section  5. 

For  the  PG  coating  used  in  the  7  inch  diameter  nozzle  insert,  the  most 
likely  modes  of  failure  are  a  compression  failure  in  the  "a”  direction  or  an 
interlaminar  shear  failure.  The  compression  strength  in  the  "a"  direction  as 
a  function  of  temperature  (Reference  6)  is  shown  in  Figure  3-23.  The  inter¬ 
laminar  shear  strength  as  a  function  of  temperature  (Reference  6)  based  on  direct 
single  shear  tests  is  shown  in  Figure  3-24.  Reference  9  gives  room  temperature 
shear  strengths  which  range  from  600  to  3200  psi  depending  on  the  method  of 
measurement.  In  evaluating  the  structural  analysis  results,  a  nominal  value  of 
1500  psi  is  used. 

For  coarse  grain  unannealed  AGSR,  ultimate  tensile  strengths  and  strains 
to  fracture  from  References  11  and  19  are  tabulated  in  Table  3-6. 
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Specific  heat  modified  to  account  for  in-depth  material  degradation 


TABLE  3-2 


MATERIAL  DENSITIES 


Material 

Density 

(Ib/ft*) 

PG  Coating 

137 

AGSR 

114 

SPEER  8882 

112 

SPEER  9139 

112 

MX4926  Carbon  Cloth  Phenolic 

89.4 

MS2600  Silica  Cloth  Phenolic 

108 

TABLE  3-3 

POISSON'S  RATIO  FOR  PYROLYTIC  GRAPHITE  COATING 
AND  THE  AGSR  GRAPHITE  SUBSTRATE 


Poisson's 

Ratio 

PG 

AGSR 

^er 

0.90 

-0.16 

0.05 

0.30 

0.07 

''zG 

— 

0.10 
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TABLE  3-4 


ESTIMATED  GRAFOIL  MECHANICAL  PROPERTIES 


Property 

Value 

Thermal  Expansion,  in/in/*F 
Across  Layers  @  70*F 

0  4000*P 

15.0  X  10“* 
70.0  X  10"‘ 

With  Layers 

S2une  as  PG 

Elastic  Modulus,  psi 

With  Layers 

0.2  X  10* 

Across  Layers 

0.4  X  10® 

Poisson's  Ratio 

Seune  as  PG 

Compressive  Strength,  psi 

12,000 

Tensile  Strength,  psi 

2,000 

TABLE  3-5 

ROOM  TEMPERATURE  YIELD  STRESSES  FOR 
AGSR  AND  PG 


Material 

Direction 

Yield  Stress 
(psi) 

PG 

"a"  (tensile) 

83,000 

"a"  (compressive) 

13,200 

"c"  (tensile) 

3,940 

"c"  (compressive) 

3,480 

Interlaminar  Shear 

300 

AGSR 

With  &  Across  Grain 

169 

Shear 

80 
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TABLE  3-6 


STRENGTH  DATA  FOR  AGSR  GRAPHITE 


Ultimate 

Strain 

Grain 

■ 

Temperature 

Tensile 

CO 

Reference 

Strength 

Fracture 

CF) 

(psi) 

(in/in) 

70 

576 

0.0015 

Across 

11 

378 

0.0014 

389 

0.0014 

490 

0.0014 

19 

800 

0.0011 

With 

19 

2000 

696 

0.0025 

Across 

11 

526 

0.0017 

482 

0.0014 

3425 

767 

0.0022 

Across 

11 

676 

0.0017 

70 

720* 

0.0028 

Radial^  . 

19 

620“ 

•  790* 

0.0018 

0.0021 

Tangential*^ 

Axial® 

^Shear  Strength 


Torsion  tests  on  radial,  tangential,  and  axial  specimens 
with  the  grain  direction  being  parallel  to  the  axis  of 
the  axial  specimen 
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SECTION  4 


THERMOCHEMICAL  ANALYSIS 

The  results  of  a  detailed  thermochemical  analysis  in  the  form  of  surface 
recession  histories  and  transient  in-depth  temperature  histories  are  required 
as  input  to  the  thermos truc^ural  analysis. 

The  thermochemical  analysis  procedure  and  the  correspondinq  design  tools 
for  generating  this  information  are  summarized  in  Section  4.1;  the  infomation 
concerning  the  nozzle  environment  which  provides  the  input  to  the  analysis  is 
presented  in  Section  4.2;  and  the  thermochemical  performance  of  the  nozzle 
materials  is  defined  in  Section  4.3. 

4.1  ■  ANALYTICAL  TECHNIQUES 

The  analysis  of  material  response  to  a  rocket  nozzle  environment  encompasses 
both  the  determination  of  surface  recession  due  to  chemical  corrosion  and  mechan¬ 
ical  erosion  and  the  evaluation  of  in-depth  heat  conduction  and  material  decompos¬ 
ition.  A  set  of  comprehensive  computerized  analysis  tools  have  been  developed 
to  define  the  thermochemical  performance  of  materials  which  are  typically 
used  in  rocket  nozzle  designs.  These  computer  codes  are  used  to  evaluate  thermal 
ablation  phenomena,  and  they  have  been  designed  to  be  applicable  for  arbitrary 
material  and  propellant  compositions  and  for  any  general  environmental  duty 
cycle  specifications. 

The  computer  codes  utilized  for  the  thermochemical  calculations  of  the 
7.0  inch  dieuneter  PG  coated  nozzle  throat  assembly  were  as  follows: 

* 

•  Aero therm  Chemical  Equilibrium  (ACE)  Program 

•  Aerotherm  Real  Gas  Energy  Integral  Boundary  Layer  (ARGEIBL)  Program 

•  Axi-Symmetric  Transient  Heating  and  Material  Ablation  (ASTHMA)  Program 

The  first  progreun  computes  the  thermochemical  state  of  any  set  of  chemical 
elemental  quantities.  For  ablative  rocket  nozzle  considerations,  the  code  assists 
in  the  computations  of  nozzle  expansion  properties,  gas  phase  transport  and  thermo¬ 
dynamic  state  properties,  pyrolysis  gas  enthalpy-temperature  relationships,  and 


The  ACE  progrcun  combines  all  the  features  of  the  Equilibrium  Surface  Thermochem¬ 
istry,  Version  3  (EST3)  code  (currently  available  to  the  rocket  technology 
community)  and  the  graphite  response  code  developed  under  Contract  F04611-C9-C- 
0081  to  account  for  surface  kinetics. 
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thermochemical  behavior  of  a  surface  material  in  contact  with  the  hyper-thermal 
boundary  layer  flow.  The  surface  removal  mechanisms  considered  by  the  code  are 
diffusion  or  kinetically  controlled  chemical  corrosion,  surface  vaporization, 
and  liquid  layer  removal. 

The  ARGEIBL  program  calculates  bulk  film  transfer  coefficients  for  non¬ 
ablating  l^^nin^u:  and/or  turbulent  boundary  layers  of  axisymmetric  bodies  or  flat 
plates  using  an  energy  integral  method.  The  solution  procedure  allows  any  real  gas 
chemical  system  because  it  accepts  generalized  thermodyneunic  property  information 
which  are  (Stained  from  ACE  code  computations.  Variations  in  axial  surface 
temperature  and  free  streeun  properties  are  accounted  for  in  the  integral  procedure. 
The  convective  heat  transfer  coefficient  distributions  computed  using  the 
ARGEIBL  technique  have  been  confirmed  by  solutions  generated  with  the  more  sophis¬ 
ticated  Boundary  Layer  Integral  Matrix  (BLIMP)  computer  code  (Reference  20) . 

Both  the  ACE  and  ARGEIBL  codes  are  utilized  to  provide  information  needed 
by  the  ASTHMA  material  response  code.  The  ASTHMA  code  is  used  to  evaluate  the 
transient  response  of  nondecomposing  materials  such  as  graphite  within  which 
two-dimensional  heat  conduction  effects  are  important.  The  surface  boundary 
condition  options  include  surface  recession  due  to  either  diffusion  or  kinetically 
controlled  chemical  corrosion.  More  complete  descriptions  of  the  capabilities 
of  the  ACE,  ARGEIBL  ,  and  ASTHMA  computer  codes  are  given  in  References  21 
through  24. 

The  sequence  of  events  in  the  analysis  procedure  is  outlined  in  Figure  4-1. 
The  major  results  of  each  analysis  step  in  the  procedure  is  defined  along  with 
how  this  result  is  used  to  obtain  the  desired  analysis  predictions.  The  analysis 
starts  by  defining  the  boundary  layer  edge  properties  from  the  flow  field  and 
ACE  expansion  analyses  and  concludes  with  the  definition  of  the  thermal 
performance  parameters,  namely,  surface  recession  and  in-depth  temperature 
distributions.  This  latter  information  is  required  as  input  to  the  structural 
analysis  which  is  discussed  in  Section  5. 

4.2  NOZZLE  ENVIRONMENT 

The  nozzle  environment  which  consists  of  the  high  speed  particle  laden 
propellant  combustion  product,  provides  the  driving  force  for  the  chemical  and/ 
or  mechani,cal  surface  recession  and  the  rise  in  the  in-depth  temperatures. 

Before  a  two-'dimensional  thermochemical  analysis  can  be  performed,  this  nozzle 
environment  must  be  quantified  in  terms  of  the  motor  ballistics  and  the  propel¬ 
lant  constituents  (Section  4.2.1),  the  boundary  layer  edge  properties  (Section 
4.2.2),  and  the  energy  transfer  mechanisms  at  the  nozzle  surface  (Section  4.2.3). 
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4.2.1  Motor  Ballistics  and  Propellant 

The  two  possible  chaunber  pressure  histories  which  were  assumed  are 
shown  in  Figure  4-2.  The  thermochemical  analysis  was  performed  for  the  chamber 
pressure  which  had  a  nominal  pressure  of  1000  psi  for  a  nominal  burn  time 
of  60  seconds.  Thermal  analysis  results  are  reported  at  motor  burn  times  of 
1.5,  10,  30,  and  60  seconds.  These  results  for  motor  burn  times  of  1.5,  30 
and  60  seconds  were  used  as  input  to  the  structural  analysis.  The  1.5  and  30 
second  results  provide  a  slightly  conservative  answer  for  the  motor  wit'i  the 
30  second  burn  time. 

The  propellant  considered  in  the  analysis  was  a  21  percent  aluminized 
PVC  composite  with  a  flame  temperature  of  5700*F  at  1000  psia.  The  constituents 
of  this  propellant  are  shown  below. 


Percent  by  Mass 


Ammonium  Perchlorate 

59.9 

Aluminum 

21.1 

Polyvinylchloride 

8.19 

DOA 

10.56 

4.2.2  Boundary  Layer  Edge  Environment 

The  boundary  layer  edge  environment  is  defined  by  performing  an  expansion 
analysis  using  the  ACE  computer  code  with  the  propellant  constituents  as  input 
data.  The  results  of  this  analysis  are  shown  in  Figures  4-3  and  4-4  as  a  function 
of  the  non-dimensional  flow  parameter  pu.  This  flow  parameter  is  the  ratio  of 
the  local  mass  flux  divided  by  the  mass  flux  at  the  sonic  point  (nozzle  throat) . 

The  boundary  layer  edge  properties  presented  in  Figures  4-3  and  4-4  are 
related  to  a  particular  nozzle  location  through  a  flow  field  analysis.  For  a 
submerged  nozzle,  the  one-dimensional  assumption  is  not  applicable  for  this 
flow  field  analysis  but  a  sophisticated  two-dimensional  analysis  was  outside  the 
scope  of  the  present  study.  In  Reference  25  a  detailed  flow  field  analysis  was 
performed  for  a  submerged  nozzle  which  was  similar  to  the  one  shown  in  Figure 
3-1.  Because  of  this  similarity,  the  flow  .'ield  results  for  the  nozzle  geometry 
presented  in  Figure  3-1  were  approximated  from  those  presented  in  Reference  25. 
These  approximate  results  are  presented  in  Figure  4-5. 

4.2.3  Wall  Surface  Energy  Balance 

For  the  problem  considered  here,  the  ASTHMA  computer  code  performs  the 
energy  balance  for  a  control  volume  fixed  to  the  receding  surface  as  indicated 
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in  the  sketch  below 
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The  energy  fluxes  leaving  the  control  volume  include  conduction  into  the  material, 
radiation  away  from  the  surface,  and  gross  blowing  at  the  surface.  Energy 
inputs  to  the  control  volume  include  radiation  from  the  combustion  products 
enthalpy  flux  due  to  the  convection  of  material  to  the  surface  associated 
with  surface  recession  and  all  the  diffusion  energy  fluxes  from  the  gas  phase 
boundary  layer.  For  a  film  coefficient  model,  the  energy  balance  becomes  for 
unequal  mass  diffusion  coefficients 


‘’e“e^H*“r  “  NrfUdge  gas  *  ^e^e^M 


(Z*  -  Z*  )h. 

le  iw  1 
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The  convection  and  mass  diffusion  coefficients  required  by  the  above  equation 
are  defined  in  Section  4. 2. 3.1,  and  the  radiation  terms  are  defined  in  Section 
4. 2. 3. 2.  The  chemical  potential  and  enthalpy  terms  are  generated  using  the 
ACE  computer  code,  and  this  information  will  not  be  presented  here.  However, 
the  thermochemical  maps  for  the  surface  materials  are  presented  in  Section  4. 2. 3. 3. 


4. 2. 3.1  Convective  Heat  Transfer  Coefficients 

The  convective  heat  transfer  coefficients  were  calculated  using  the 
Aerotherm  Real  Gas  Energy  Integral  Boundary  Layer  (ARGEIBL)  program.  This  pro¬ 
gram  has  been  used  previously  to  perform  boundary  layer  analysis  of  the  type 
presented  here  (i.e..  References  5,  20,  and  25).  The  primary  inputs  for  this 
program  are 

e  The  real  gas  (With  the  Al20^  particles  removed)  Mollier  tables  for 
the  propellant. 

e  The  boundary  layer  edge  properties  as  a  function  of  pu. 

e  The  pu  ratio  as  a  function  of  nozzle  location. 

The  first  two  items  of  input  data  were  generated  using  the  Aerotherm  Chemical 
Equilibrium  (ACE)  program  and  the  third  item  was  obtained  from  Figure  4-5. 


4-4 


One  other  item  to  note  is  that  the  initiation  of  the  boundary  layer  (X  =  0)  is 
defined  in  Figure  3-1  and  is  on  the  backface  of  the  submerged  nozzle.  The 
convective  heat  transfer  coefficient  results  from  the  ARGEIBL  prograun  are 
presented  in  Figure  4-6  for  two  wall  temperatures.  Based  on  previous  analyses 
(References  5,  20,  and  25)  these  coefficients  were  multiplied  by  0.75  prior  to 
using  them  in  the  two-dimensional  heat  conduction  analysis.  The  mass  transfer 
coefficients  were  generated  from  these  convective  heat  transfer  coefficients 
using  the  relation 

where  P^u^Cw  -  mass  transfer  coefficient 
e  e  M 

^e'^e^H  ”  heat  transfer  coefficient 
Le  -  Lewis  number 

4. 2. 3. 2  Radiation  Heating 

This  section  is  concerned  with  the  two  radiation  terms  in  liquation  (4-1)  . 
The  radiant  flux  absorbed  by  the  surface  (4j.ad  in^  defined  as  the  incident 
radiation  multiplied  by  the  absorptance  of  the  wall.  The  radiant  flux  emitted 
by  the  surface  ^y^)  is  defined  as  that  emitted  by  a  black  body  multiplied 

by  an  emittance  factor  and  by  a  view  factor.  Therefore  the  net  radiation  to 
the  surface  is  defined  by 

^net  rad  w^inc  rad  w  w 

where 

o  ,e  -  absorptivity  and  emissivity,  respectively,  of  wall  material 
w  w 

(a^  and  are  equal  by  Kirchoff's  Law  for  a  gray  body) 
o  -  Stefan-Boltzmann  constant 
T^  -  wall  temperature 
F  -  view  factor 

In  order  to  model  the  radiation  transport  between  a  particle  laden  streeun 
and  the  nozzle  surface,  it  was  assumed  that  the  stream  was  optically  thick  and 
that  the  particles  and  wall  exchange  radiant  energy  as  if  they  were  two  parallel 
plates.  In  this  way  multiple  reflections  between  the  wall  and  stream  were  taken 
into  account.  In  addition,  it  was  assumed  that  both  the  stream  and  wall  behave 
as  gray  bodies  and  that  they  emit  and  reflect  radiant  energy  diffusely.  Based 
on  the  above  assumption,  the  net  radiant  heat  flux  relation  is  given  as 
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(4-4) 


where 


%et  rad  *  =e£f  <'”'a  "  ‘”'w> 

(effective  emissivity)  = 

Eg  -  particle  laden  stream  emissivity 
Tg  -  local  stream  static  temperature 


(4-5) 


To  define  the  net  radiation  heat  flux  using  Equation  (4-4)  the  local 
static  stream  temperature  and  the  effective  emissivity  had  to  be  evaluated.  The 
stream  temperature  was  obtained  from  the  isentropic  expansion  data  presented 
previously  in  Figures  4-3  and  4-4.  The  effective  emissivity  was  evaluated  from 
Equation  (4-5)  using  a  value  of  0.9  for  the  wall  emissivity  and  using  the  follow¬ 
ing  relationship  (Reference  5)  for  the  emissivity  of  the  particle  laden  combustion 
products . 


(4-6) 


where 

C  -  experimentally  determined  constant  (0.808  for  the  propellant  being 
considered  in  this  analysis) 

n  -  percentage  of  aluminum  loading 

p  -  local  density  of  propellant  combustion  species  (Ib/ft’) 

D  -  local  beam  length,  usually  taken  as  the  diameter  (in.) 


4.2. 3. 3  Surface  Thermochemistry  Maps 

In  order  to  perform  a  transient  two-dimensional  thermal  analysis  of  the 
nozzle  throat  assembly,  surface  thermochemistry  maps  are  required  for  the  Speer 
8882  and  Speer  9139  bulk  graphites,  the  PG  coating,  and  the  MX  4926  carbon  cloth 
phenolic.  The  surface  thermochemical  response  of  the  first  three  materials  is 
controlled  by  the  finite  (as  opposed  to  equilibrium)  reaction  rates  of  the  sur¬ 
face  carbon  atoms  with  the  H2O,  CO2,  and  K2  components  of  the  solid  propellant. 

The  kinetically  controlled  reaction  rate  models  developed  under  Air  Force  Contract 
F04611-69-C-0081  (Reference  24)  were  utilized  to  compute  the  surface  response 
of  these  graphite  materials,  and  solutions  were  performed  using  the  ACE  code. 
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The  kinetic  response  model  incorporates  a  Langmuir  type  surface  reaction 
inhibition  expression  into  an  Arrenhius  temperature  dependent  rate  relation. 

The  rate  expression  and  reaction  kinetic  rate  constants  are  given  in  Tedsle  4-1. 

The  partial  pressure  quantities  were  evaluated  in  the  ACE  code  using  the  film 
coefficient,  multicomponent  unequal  diffusion  boundary  layer  model.  The  area 
ratio,  local  static  pressure  and  temperature,  and  the  local  heat  and  mass  transfer 
coefficient  value  for  which  the  kinetically  controlled  response  calculations 
were  performed  are  given  in  Table  4-2.  The  non-dimensional  kinetically  controlled 
surface  thermochemistry  maps  are  shown  in  Figure  4-7.  For  comparison,  the 
non-dimensional  diffusion  controlled  map  is  also  shown  in  this  figure. 

The  steady  state  surface  thermochemistry  map  for  the  MX4926  carbon  cloth 
phenolic  is  shown  in  Figure  4-8.  For  the  particular  environment  which  it  is 
subjected  to,  the  steady  state  assumption  is  probably  reasonable  for  this  material. 
However,  the  primary  reason  for  this  assumption  is  that  the  two-dimensional 
computer  code  does  not  account  for  material  degradation,  and  thus,  it  can  not 
account  for  variable  pyrolysis  gas  rates. 

4.3  NOZZLE  WALL  MATERIAL  THERMOCHEMICAL  PERFORMANCE 

A  two-dimensional  thermal  analysis  including  the  effects  of  surface  abla¬ 
tion  was  performed  using  the  Axi -Symmetric  Transient  Heating  and  Material  Abla¬ 
tion  (ASTHMA)  program  (Reference  21)  .  An  analysis  was  performed  using  two  nodal 
layouts:  the  nozzle  insert  nodal  grid  (Figure  4-9)  and  the  nozzle  assembly  nodal 

grid  (Figure  4-10)  .  These  grids  represent  the  two  possible  extremes  with  1)  the 
validity  of  the  nozzle  insert  nodal  grid  dependent  upon  the  assumption  of  infinite 
thermal  resistance  between  the  insert  and  the  adjacent  materials,  and  2)  the 
validity  of  the  nozzle  assembly  grid  dependent  upon  the  assumption  of  zero 
thermal  resistance  between  the  insert  and  the  adjacent  materials.  Because  of 
the  fabrication  procedure,  tolerances,  and  bonding  materials,  the  latter  assump¬ 
tion  is  believed  to  be  more  realistic. 

The  surface  temperatures,  surface  recession  rates,  and  net  surface 
recessions  are  presented  in  Figures  4-11  through  4-13  for  the  nozzle  insert  grid 
and  for  the  entire  60  second  motor  firing  and  in  Figures  4-14  through  4-16  for 
the  nozzle  assembly  grid  and  for  the  first  30  seconds  of  the  motor  firing.  The 
results  presented  in  these  figures  were  based  on  the  original  "c"  direction 
thermal  conductivity  values  for  pyrolytic  graphite.  The  PG  results  presented 
in  Figure  4-11  through  4-16  which  are  compared^le  show  that  the  predicted  perform¬ 
ance  of  the  PG  coating  is  independent  of  the  nodal  grid. 

The  predicted  performance  presented  in  Figures  4-14  through  4-16  for  both  the 
downstream  bulk  graphite  and  the  MX4926  carbon  cloth  phenolic  materials  indi¬ 
cate  no  design  problem.  However,  the  predicted  performance  of  the  upstream 
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Speer  8882  bulk  graphite  material  is  marginal  for  a  60  second  firing.  The 
problem  is  that  the  predicted  surface  recession  for  this  material  would  be 
approximately  320  mils  for  a  60  second  firing.  This  amount  is  larger  than 
the  combined  thickness  of  the  Speer  8882  lip  (200  mils)  and  the  PG  coating 
(60  mils) .  A  comparison  of  this  thickness  with  the  predicted  recession  meems 
that  the  interface  between  the  PG  coating  and  the  AGSR  substrate  would  be  exposed. 
Because  this  interface  would  be  partially  shielded  as  shown  by  the  sketch  below, 
the  actual  recession  at  the  interface  would  be  less  than  the  predicted  value  of 
320  mils.  However,  this  shielding  may  not  be  sufficient  to  prevent  undercut¬ 
ting  of  the  PG  coating. 
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Isotherm  temperature  contours  for  motor  firing  times  of  1.5,  5.0,  10.0, 
30.0,  and  60.0  and  for  the  nozzle  insert  grid  are  presented  in  Figures  4-17 
through  4-21.  These  figures  show  the  gradual  progression  of  the  isotherms  into 
the  AGSR  as  the  motor  firing  continues.  For  example,  at  t  >  1.5  seconds,  the 
interface  of  the  PG  coating  and  the  AGSR  substrate  is  approximately  1000 "R. 

At  t  «  1C  seconds,  this  interface  has  increased  to  2000®R  and  the  majority  of 
the  substrate  is  above  1000‘R.  At  t  «  30  seconds,  the  interface  is  approximately 
3000*R,  and  at  t  =  60  seconds,  it  is  approximately  40o0°R.  In  addition,  these 
figures  show  that  the  upstreeun  portion  of  the  substrate  is  at  higher  temperatures 
than  either  the  throat  or  downstream  sections.  This  is  due  to  the  substrate 
being  relatively  thinner  at  the  upstream  section. 

Isotherm  temperature  contours  for  motor  firing  times  of  1.5,  5.0,  10.0 
and  30.0  seconds  and  for  the  nozzle  assembly  grid  are  presented  in  Figures  4-22 
through  4-25.  For  the  nozzle  insert,  the  primary  difference  between  the  tempera¬ 
ture  contours  presented  in  these  figures  and  the  corresponding  ones  presented 
previously  is  the  progression  towards  the  throat  of  the  upstrecun  and  downstream 
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portions  of  the  isotherms.  This  progression  is  caused  by  thermal  energy  being 
transmitted  through  the  upstream  and  downstream  bulk  graphites  and  into  the 
AGSR  substrate  and  is  a  result  of  the  thermal  conductivity  of  the  bulk  graphites 
being  much  higher  than  for  the  "c"  direction  of  the  PG  coating.  The  differences 
in  the  thermal  profiles  at  the  entrance,  throat,  and  exit  of  the  nozzle  insert 
for  the  two  nodal  grids  are  illustrated  in  Figure  4-26.  A  comparison  of  these 
profiles  shows  the  following: 

•  For  the  entrance  and  exit  sections  of  the  AGSR  the  predicted  tempera¬ 
tures  for  t  =  30  seconds  using  the  nozzle  assembly  grid  are  approxi¬ 
mately  the  same  as  those  predicted  for  t  =  60  seconds  using  the 
nozzle  insert  grid 

•  At  30  seconds,  the  predicted  temperatures  at  the  throat  section  are 
approximately  the  same  for  both  nodal  grids. 

The  differences  between  the  nozzle  insert  temperatures  for  the  two  nodal  grids 
do  not  significantly  affect  the  thermal  performance  of  the  throat.  However, 
these  temperature  differences  might  have  a  significant  effect  on  the  structural 
performance  of  the  insert,  and  this  effect  will  be  investigated  in  the  structural 
analysis  which  is  presented  in  Section  5. 

A  two-dimensional  ana''ysis  was  performed  using  the  revised  "c"  direction 
PG  thermal  conductivity  properties  and  the  nozzle  insert  nodal  grid.  By  comparing 
the  results  of  this  analysis,  which  are  presented  in  Figures  4-27  through  4-31, 
with  those  presented  in  Figures  4-17  through  4-21,  one  notes  that  the  revised 
PG  properties  cause  a  marked  increase  in  the  thermal  penetration  of  the  AGSR  sub¬ 
strate.  This  increase  is  best  seen  by  comparing  the  location  of  the  3000‘’R 
and  4000‘’R  isotherms  presented  in  Figures  4-21  and  4-31.  The  primary  reason  for 
performing  the  analysis  using  the  revised  PG  properties  was  to  show  the  effect 
that  the  "c"  direction  PG  thermal  conductivity  has  on  the  AGSR  substrate  tempera¬ 
tures  . 
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Table  4-1 

KINETICALLY  CCNTROLl.KD  REACTION  RATE  EXPRESSIONS  USED  FOR  PG  lAYER 
AflD  SPEER  BULK  GRAPHITE  SURFACE  THERMOCHEMISTRY  SOLUTIONS* 


Material 

Pre-Exponent  Factors 
(units  to  give  X^,,  below) 

Activation  Energies 
(cal/mole) 

®2 

=1 

=2 

Layer  PG 

12.5 

0.77 

46,000 

55,500 

Edge  PG*» 

1.65  X  10^ 

4.12  X  10* 

65,500 

129,500 

ATJ  (BulkJ  Graphite 

1,98  y.  10* 

4.94  X  10* 

i 

1 

_ _ Biexp(-E^/RT^)  [■/  _  ^COW  (  _  P^J. 

C*  (1.  +  +  PcOj  *  Pco^  ^  "2°  *^1120  j  ^^02/ 


moles  oZ  reaction/ft' cec 


b 


No  edge  PG  solutions  wore  performed;  constants  included  for  completeness. 


SUMMARY  OF  EDGE  CONDITIONS  FOR  SURFACE  THERMOCHEMISTRY  MAPS 
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FIGURE  4-1.  FLOW  CHART  FOR  AEROTHERMOOYNAMIC  ANALYSIS 
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Figure  4-10.  Nozzle  Assembly  Nodal  Grid 
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igure  4-17.  Isotherni  Ter.perature  Profiles  for  Nozzle  Insert 
Nodal  Grid  -  t  =  1.5  seconds 


i  I 


^  V)  V)  >  I\j  \ 

>|c  X+  <3  O  □ 


4  i 

!  i 
I  : 

I  I  / 


'  I  i 


/  /■/' 


'  f 

fi  4 

'  :,:i 


4-29 


Figure  4-18.  Isotherm  Temperature  Profiles  for  Nozzle  Insert 
Nodal  Grid  -  t  =  5.0  seconds 
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Figure  4-21.  Isotherm  Temperature  Profiles  for  Nozzle  Insert 
Nodal  Grid  -  t  =  60.0 
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SECTION  5 


STRUCTURAL  ANALYSIS 

Thtt  overall  objective  of  the  structural  analysis  was  to  define  the  struc¬ 
tural  adequacy  of  the  large  PG  coated  nozzle  throat  insert  using  available  de¬ 
sign  techniques  and  material  properties.  To  evaluate  this  objective  required 
performing  structural  analyses  at  three  stages  during  the  fabrication  of  the  PG 
coated  insert  and  at  several  times  during  the  motor  firing.  As  outlined  in 
Table  5-1,  structural  analyses  were  performed  for 

e  The  heat  up  of  the  AGSR  substrate  to  2000*C 

e  The  cooldoim  of  the  PG  coated  insert  from  2000*C  to  room  temperature 

e  The  machining  of  the  insert 

e  The  motor  firing  corresponding  to  times  of  1.5,  30.0,  and  60.0 
seconds 

In  addition,  the  results  from  the  structural  analyses  were  to  provide  some 
guidelines  concerning  the  Influence  of  the  following  design  and/or  analysis 
variables. 

e  Define  the  effect  on  the  stress  (or  strain)  distribution  resulting 
frOTi  the  placement  of  grooves  in  the  upstreeun  and  downstream  ends 
of  the  AGSR  substrate 

e  Define  the  effect  on  the  predicted  stress  (or  strain)  distribution 
tions  resulting  from  the  Inclusion  of  plasticity  effects  in  the 
analysis 

e  Define  the  effect  on  the  predicted  stress  (or  strain)  distributions 
resulting  from  using  the  nozzle  insert  or  nozzle  assembly  tempera¬ 
tures 

e  Define  the  effect  on  the  predicted  stress  (or  strain)  distributions 
resulting  from  variations  in  the  backside  boundary  constraints. 

As  shown  by  Table  5-1,  the  effect  of  the  grooves  in  the  substrate  were  evalu¬ 
ated  for  the  cooldown  of  the  PG  coated  insert  and  for  a  motor  firing  time  of 
1.5  seconds,  the  effect  of  including  plasticity  in  the  analysis  was  evaluated 
only  for  the  cooldown  of  the  Insert,  the  effect  of  the  temperature  distribution 
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was  evaluated  for  a  motor  firing  time  of  30  seconds,  and  the  effect  of  the  back¬ 
side  constraint  was  evaluated  for  motor  firing  times  of  1.5,  30.0,  and  60.0 
seconds. 

The  material  properties  used  in  the  structural  analyses  have  been  pre¬ 
sented  in  Section  3.2.2.  The  analytical  tool  used  for  these  structural  analy¬ 
ses  was  the  Deformation  Orthotropic  Axlsymmetric  Solution  of  Inelastic  Solids 
(DOASIS)  computer  progr€un,  and  it  is  described  in  Section  5.1.  The  geometric 
representation  of  the  nozzle  insert  is  outlined  in  Section  5.2,  and  the  analyti¬ 
cal  modeling  of  the  insert  fabrication  process  is  presented  in  Section  5.3.  The 
insert  external  boundary  conditions  are  defined  in  Section  5.4,  and  the  in-depth 
thermal  loading  is  outlined  in  Section  5.5.  Finally,  the  thermostructural  re¬ 
sponse  of  the  insert  for  the  fabrication  cooldown  and  motor  firing  environments 
is  described  in  Section  5.6. 

5.1  ANALYTICAL  TECHNIQUES 

The  primary  analytical  tool  used  in  the  structural  analysis  of  the  rocket 
nozzle  insert  was  the  DOASIS  computer  code.  The  DOASIS  program  is  a  very  gen¬ 
eral  computer  code  which  has  the  capability  of  analyzing  the  elastic-plastic 
behavior  of  either  axlsymmetric  bodies  of  revolution,  plane  stress,  or  plane 
strain  solids.  The  program  is  capable  of  handling  generalized  multi-modulus, 
orthotropic  materials.  The  influences  of  thermal,  pressure,  shear,  concentrated 
force,  centripetal  force,  body  force,  and  pre-stress  loadings  ere  accounted  for. 
In  addition,  displacement  and  interference  fit  boundary  conditions  are  available 
options . 

DOASIS  obtains  the  above-mentioned  solutions  by  use  of  the  finite  element 
method.  That  is,  a  continuous  body  is  subdivided  into  small,  discrete,  quadri¬ 
lateral  or  triangular  elements.  The  peripheral  program  MESHGEN  generates  this 
assemblage  of  elements  for  DOASIS  by  a  Laplacian  interpolation  technique.  The 
DOASIS  program  then  uses  this  mesh  and  the  material  properties,  loading  and 
boundary  condition  data  to  generate  a  stiffness  matrix  for  the  given  program. 

The  elastic-plastic  problem  is  then  solved  by  the  most  efficient  method  avail¬ 
able,  i.e.,  by  using  an  iterative  technique  in  which  only  the  nodal  point  force 
vector  changes  during  the  iterative  process.  Hence  the  stiffness  mcitrix  of  the 
structure  need  only  be  generated  once  and  decomposed  once. 

The  program  allows  a  general  treatment  of  orthotropy.  The  material  prop¬ 
erties  in  each  of  the  principal  directions  (designated  as  the  1-2-6  axis  system) 
are  independent  and  are  input  in  tabular  form  as  arbitrary  functions  of  tempera¬ 
ture.  The  sketch  below  provides  a  convenient  representation  of  the  generality. 
The  program  allows  the  principal  (1-2)  axes  to  be  inclined  at  an  angle  (fi  to  the 
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r-z  axis  system.  The  angle  (P  can  vary  from  point  to  point  in  an  arbitrary  man¬ 
ner.  The  meridianal  (r-z  plane)  geometry  can  be  completely  arbitrary ,  but  the 
loading  must  be  axisymmetric . 


For  plasticity  problems,  the  code  uses  the  associated  flow  rule  of  plas¬ 
ticity  theory.  The  magnitude  of  any  plastic  strain  increment  is  determined  by 
the  strain  hardening  portion  of  the  stress-strain  curves  (the  tangent  moduli) 
and  their  directions  (ratio  of  the  components  with  respect  to  one  another)  by 
the  outward  normal  to  the  yield  surface  as  described  by  the  current  stresses. 
The  yield  surface  used  here  is  a  modified  form  of  the  Hill-Hu  orthotropic  yield 
surface  given  by 
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where  "a  =  equivalent  stress,  Y  »  equivalent  yield  stress  and  -  anisotropic 
pareuneters.  A  full  development  of  the  theory  can  be  found  in  the  Aerotherm 
DOASIS  User's  Manual  (Reference  15). 
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Th«  peripheral  programs  besides  HESHGEN  which  feed  input  data  to  DOASZS 
are  (1)  TEMINTi  a  generalized  temperature  interpolator  program  idiich  takes  an 
arbitrary  two-dimensional  temperature  distribution  and  interpolates  to  find  the 
temperatures  of  all  the  points  required  by  the  DOASIS  code;  (2)  GETPR,  a  progreun 
designed  to  translate  a  given  pressure  distribution  on  the  outside  surface  of  a 
sphere-cone  geometry  to  the  outside  surface  of  a  similar  structural  finite  ele¬ 
ment  mesh  and  generates  a  general  "cosine"  distribution  of  pressure  for  the 
backside  surface  of  the  structural  finite  element  mesh;  (3)  INTERG,  a  program 
which  uses  the  results  of  a  DOASIS  computer  run  to  calculate  maximum  and  minimum 
stresses,  strains  and  strain  ratios,  and  determines  the  locations  of  all  the  com¬ 
ponents  of  the  different  variables;  and  (4)  CONTOUR,  a  program  which  is  designed 
to  make  contour  plots  of  the  analysis  results  generated  by  the  DOASIS  computer 
code. 


5.2  GEOMETRIC  REPRESENTATION  OF  ROCKET  NOZZLE  INSERT 

5.2.1  Insert  Without  Stress  Relieving  Grooves 

The  finite  element  structural  representation  of  the  rocket  nozzle  insert 
without  stress  relieving  grooves  is  shown  in  Figure  5-1.  A  total  of  195  ele¬ 
ments  were  used  to  represent  the  AGSR  portion  of  the  insert,  and  45  elements 
accounted  for  the  PG  coating.  A  finer  mesh  was  used  in  the  PG  coating  because 
of  the  large  stress  and  strain  gradients  expected  there.  The  dimensions  of  the 
entire  mesh  correspond  to  the  geometry  of  the  insert  before  machining  0.1  inch 
from  the  ends.  The  structural  mesh  of  the  insert  without  grooves  was  used  pri¬ 
marily  to  compute  cooldown  stresses.  These  computations  were  used  as  a  basis 
for  determining  the  effect  of  grooves  in  reducing  the  residual  stresses 
in  an  insert  resulting  from  cooldown.  Only  one  motor  firing  structural  analy¬ 
sis  was  performed  using  the  mesh  without  grooves.  In  that  particular  case,  the 
mesh  was  modified  by  shortening  the  end  elements  to  account  for  the  machining 
process  which  was  performed  prior  to  placing  the  insert  in  the  rocket  nozzle 
structural  shell. 

5.2.2  Insert  with  Stress  Relieving  Grooves 

The  finite  element  structural  representation  of  the  rocket  nozzle  insert 
with  stress  relieving  grooves  is  shown  in  Figure  5-2.  A  total  of  391  elements 
were  used  to  represent  the  AGSR  portion  of  the  insert,  and  45  elements  accounted 
for  the  PG  coating.  A  finer  mesh  was  used  with  this  geometry  for  several  rea- 
sonf : 

1.  The  presence  of  the  grooves  results  in  larger  stress  concentration 
gradients  near  the  ends  of  the  insert 
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2.  El«m«nts  at  tha  azaSs  %f«ra  raqulrad  to  ba  0.1  inch  thick  so  that  thay 
could  ba  raaovad  tdian  analyses  wars  parfomad  on  a  awchinad  insert. 

3.  Motor  firing  loadings  and  boundary  conditions  %iara  applied  to  tha 
insert  with  grooves,  therefore  large  radial  thersMl  gradients  ware 
present  in  the  PG  coating. 

4.  Tha  P6  coating  ablated  during  tha  motor  firing,  hence  a  variable 
number  of  elements  were  removed  from  the  coating  depending  on  the 
number  of  seconds  elapsed  in  the  motor  firing. 

5.3  MODELING  OF  THE  INSERT  FABRICATION  PROCESS 

The  fabrication  of  the  rocket  nozzle  insert  was  modeled  by  decomposing 
the  process  into  several  steps.  Each  step  was  chosen  such  that  there  was  no 
need  for  any  simplifying  assumptions  with  regard  to  the  material  response.  The 
results  of  each  preceding  step  were  inserted  as  input  data  to  the  next  step  in 
the  process,  if  necessary.  The  fabrication  process  was  segmented  as  follows 
(also  refer  to  Figure  5-3) : 

1.  An  unmachined  billet  of  AGSR  Without  a  PG  coating  and  in  a  stress 
free  state  is  machined  so  that  the  stress  relieving  grooves  are 
placed  on  both  ends  of  the  substrate,  and  the  nozzle  contour  is 
placed  on  the  inner  surface  of  the  substrate. 

2.  The  billet  of  AGSR  is  next  heated  from  70°  to  2000*C. 

3.  The  PG  coating  is  placed  on  the  inner  surface  of  the  AGSR  substrate 
at  3634»F  (2000*0. 

4.  The  combined  PG  coating  and  AGSR  substrate  are  then  cooled  to  70*F. 

5.  The  Insert  is  then  machined  so  that  0.1  inch  is  removed  from  each 
end. 

This  completes  the  fabrication  process.  In  each  step,  the  surface  of  the 
insert  was  assumed  to  be  stress  free. 

5.4  DEFINITION  OF  INSERT  BOUNDARY  CONDITIONS  FOR  MOTOR  FIRING 

Once  the  fabrication  process  is  complete,  the  insert  is  placed  in  posi¬ 
tion  in  the  rocket  nozzle  structural  shell.  At  this  point  it  becomes  very  im¬ 
portant  to  properly  define  the  boundary  conditions  on  the  exterior  surfaces  of 
the  insert.  These  boundary  conditions  are  illustrated  in  Figure  5-4.  The  vari¬ 
ous  faces  of  the  insert  are  defined  as  the  inner,  outer,  upstream,  or  downstream 
surfaces.  Each  surface  will  be  considered  separately. 
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5.4.1  lnn«r  Surface  Boundary  Condition 

Th«  innnr  surfac*  of  th«  insert  faces  the  nozzle  chasibur  environment.  It 
is,  therefore,  exposed  to  the  severe  pressure,  thermal,  and  corrosive  environ* 
ment  of  the  propellant  combustion  products.  For  the  7-inch  diasMter  PG  coated 
nozzle  insert,  the  pressure  varies  approximately  linearly  from  57.5  atznspheres 
at  the  upstream  edge  to  26.7  atmospheres  at  the  downstream  edge.  The  severe 
thermochemical  environment  causes  the  PG  coating  to  ablate  so  that  the  thickness 
of  the  coating  varies  during  the  motor  firing. 

5.4.2  Outer  Surface  Boundary  Condition 

The  outer  surface  of  the  AGSR  is  overwrapped  with  0.006  inches  of  teflon 
which  is  separated  from  the  silica  cloth  phenolic  by  a  0.10  inch  thick  bondline 
of  EA913  adhesive.  Various  boundary  conditions  were  placed  on  this  surface  de¬ 
pending  on  the  particular  case  considered.  In  cases  where  the  1.5  second  sK>tor 
firing  conditions  tfere  considered,  it  was  assumed  that  the  bondline  and  teflon 
%iere  rigid  compared  with  the  AGSR,  and  hence  there  was  do  displacement  in  the 
radial  direction  but  free  movement  in  the  axial  direction.  In  cases  where  the 
30  and  60  second  motor  firing  conditions  were  considered,  it  was  assumed  that 
the  teflon  evaporated  in  one  care,  and  that  both  the  teflon  and  the  bondline 
evaporated  in  another  case.  Therefore,  the  boundary  condition  on  the  outer  sur¬ 
face  in  these  two  cases  was  a  stress  free  condition  until  the  insert  expanded 
to  contact  the  silica  cloth  phenolic.  This  is  also  referred  to  as  an  interfer¬ 
ence  fit  boundary  condition.  In  another  case,  a  0.2  inch  thickness  of  multi¬ 
layered  grafoil  was  assumed  to  exist  between  the  AGSR  and  the  silica  cloth 
phenolic.  In  this  case,  the  outer  surface  of  the  grafoil  was  assumed  to  be  a 
rigid  boundary  but  free  axial  movement  was  allowed.  The  inner  surface  of  the 
grafoil  was  assumed  bonded  to  the  outer  surface  of  the  AGSR. 

5.4.3  Upstream  Surface  Boundary  Condition 

The  upstream  surface  of  the  rocket  nozzle  insert  is  separated  from  the 
Speer  8882  material  by  a  0.005  to  0.02  inch  thick  bondline  consisting  of  Arm¬ 
strong  A-2.  This  bondlir.e  was  assumed  to  be  structurally  weak,  hence  the  pres¬ 
sure  at  the  inner  edge  of  this  surface  was  assumed  to  be  transmitted  throughout 
the  Armstrong  A-2.  Therefore,  the  boundary  condition  on  the  upstream  surface 
was  always  assumed  to  be  a  pressure  of  57.5  atmospheres. 

5.4.4  Downstream  Surface  Boundary  Condition 

The  boundary  condition  on  this  surface  was  the  most  complicated  because: 
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1.  TIm  pc«taur«  on  th«  lurfac*  must  be  greater  than  the  local  gas  static 
pressure  in  order  to  place  the  insert  in  force  equilibrium  during  the 
■otor  firing. 

2.  The  differential  surface  recession  rate  between  the  PG  coating  and  the 
bulk  graphite  continually  changes  the  area  over  which  this  reaction 
pressure  is  applied. 

3.  The  0.005  to  0.020  inch  thick  bondline  of  Armstrong  A-2  between  the 
AG8R  and  the  Speer  9139  gradually  degrades  during  the  motor  firing, 
and,  thus,  at  any  instant  of  time,  it  has  a  variable  stiffness  in  the 
radial  direction. 

After  1.5  seconds  into  the  motor  firing,  the  reaction  pressure  which  acts  over 
the  entire  surface  was  calculated  to  be  57.6  atmospheres.  However,  for  the  30 
and  60  second  cases,  the  possibility  exists  that  the  combined  effect  of  items 
2  and  3  could  reduce  the  effective  area  over  which  the  reaction  pressure  is  ap¬ 
plied.  In  this  event,  the  effective  area  is  assumed  to  include  the  downstre2un 
surface  area  between  the  inner  edge  of  the  groove  and  the  outside  radius  of  the 
substrate  (refer  Figure  5-4).  For  this  effective  surface  area,  the  average  pres¬ 
sure  is  based  on  the  local  static  pressure  of  26.7  atmospheres  being  applied  to 
the  remainder  of  the  downstream  surface  which  includes  the  area  from  the  inner 
edge  of  the  PG  coating  to  the  inner  edge  of  the  groove.  However,  because  the 
bondline  has  significantly  degraded  at  the  later  times,  the  AGSR  is  probably  in 
contact  with  the  Speer  9139  bulk  graphite.  Hence,  the  boundary  condition  over 
the  effective  area  portion  of  the  downstream  surface  for  the  30  and  60  second 
cases  was  chosen  as  one  with  no  axial  displacement  but  free  radial  movement. 

The  other  boundary  condition  used  for  the  30  and  60  second  cases  was  the  same  as 
used  for  the  1.5  sec  cases. 

5.5  DEFINITION  OF  INSERT  INTERNAL  LOADING  CONDITIONS 

The  only  Internal  loading  on  the  rocket  nozzle  insert  during  a  motor  fir¬ 
ing  is  that  resulting  from  the  two-dimensional  temperature  distribution.  The 
temperature  distribution  in  the  insert  as  a  function  of  time  was  determined  by 
the  ASTHMA  code.  However,  as  discussed  in  Section  4.3,  two  different  tempera¬ 
ture  distributions  were  calculated  depending  on  whether  the  nozzle  insert  alone 
was  investigated  or  whether  the  entire  nozzle  assembly  was  included  in  the  heat 
transfer  analysis.  For  roost  cases,  the  nozzle  insert  temperature  distributloi 
was  used  for  the  structural  analyses;  however,  a  check  case  at  30  seconds  was 
performed  with  the  entire  nozzle  assembly  temperature  distribution. 


5.6  PREDICTED  STRUCTURAL  RESPONSE 

Based  on  the  previously  described  geometry,  material  properties,  loading, 
and  boundary  conditions,  structural  analyses  were  performed  to  determine  the 
stress-strain-displacement  state  of  the  rocket  nozzle  insert  during  both  the 
fabrication  process  and  a  60  second  motor  firing.  As  outlined  in  detail  in  Table 
5-2,  calculations  were  performed  for  a  total  of  18  different  conditions.  In 
this  table,  these  conditions  are  defined  by  case  number,  solution  type  (elastic- 
plastic)  ,  the  AGSR  substrate  type  (with  or  without  stress  relieving  grooves) , 
the  temperature  of  the  Insert,  the  time  elapsed  since  the  start  of  the  motor 
firing  (unless  the  condition  was  related  to  some  state  during  the  fabrication 
process  of  the  insert) ,  the  surface  boundary  conditions  which  are  explained  in 
more  detail  in  Table  5-3,  the  reference  stress  free  state,  and  any  necessary 
explanatory  remarks.  The  results  of  these  calculations  are  subdivided  into  two 
categories:  those  relevant  to  the  fabrication  process  which  are  presented  in 
Section  5.6.1  and  those  relevant  to  the  motor  firing  which  are  presented  in  Sec¬ 
tion  5.6.2. 

5.6.1  Fabrication  Process 

As  discussed  in  Section  5.3,  the  fabrication  of  the  rocket  nozzle  insert 
can  be  subdivided  into  several  steps.  The  results  of  the  structural  analyses 
for  each  of  these  steps  is  presented  below. 

5. 6. 1.1  Heating  to  3634«F 

The  first  step  consisted  of  heating  a  billet  of  AGSR  from  70"F  to  3634 ®F 
Results  of  the  numerical  calculations  (Cases  1  and  2)  indicated  that  the  stresses 
and  strains  at  3634 ®F  were  negligible  although  the  insert  expanded  22  mils  radi¬ 
ally  outward  and  14.5  mils  axially. 

5. 6. 1.2  Cooling  to  70*F  with  PG  added 

The  second  step  in  the  fabrication  process,  from  the  structural  point  of 
view,  consisted  of  coating  the  AGSR  substrate  at  3e4G°F  with  the  pyrolytic 
graphite  coating  then  cooling  the  combination  briok  to  70®F.  Due  to  the  dif¬ 
ferences  in  the  coefficients  of  thermal  expansion  of  the  two  materials,  the  cool¬ 
ing  of  the  coated  substrate  from  3646®F  to  70®F  generates  internal  stresses. 

Cases  3,  5,  and  6  of  the  numerical  calculations  define  this  stress  state  under 
three  different  conditions.  Cases  3  and  5  assume  the  material  behavior  can  be 
characterized  as  linear  elastic,  but  Case  3  is  for  a  grooved  substrate  and  Case 
5  is  for  a  non-grooved  substrate.  Case  6  assumes  that  the  insert  has  stress  re¬ 
lieving  grooves,  but  that  the  material  behavior  can  be  characterized  as  elastic- 
plastic.  Each  of  these  cases  are  discussed  in  more  detail  below. 
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5. 6. 1.2.1  Insert  with  Grooves  vs.  Insert  without  Grooves 

The  basic  reason  for  the  use  of  stress  relieving  grooves  in  the  insert 
'  was  to  reduce  the  large  interlaminar  shear  stresses  in  the  insert  at  the  inter¬ 
face  between  the  AGSR  and  the  PG  coating.  Figure  5-5  illustrates  to  what  ex¬ 
tent  these  shear  stresses  were  reduced  by  the  grooves.  The  results  indicate 
that  the  grooves  reduce  the  peak  shear  stresses  on  the  upstream  side  of  the 
insert  by  30  percent,  and  on  the  downstream  side  by  about  27  percent.  Further 
results  (Figures  5-6  and  5-7)  indicate  that  the  hoop  stresses  in  the  PG  coating 
for  the  grooved  substrate  are  reduced  by  as  much  as  21  percent.  However,  for 
^  certain  regions  in  the  AGSR,  the  hoop  stresses  and  strains  are  larger  in  the  in¬ 
sert  with  grooves.  The  results  are  illustrated  in  Figures  5-8,  5-9,  and  5-10. 
Furthermore,.  Figures  5-11  and  5-12  indicate  regions  in  the  AGSR  where  failure  of 
the  material  would  be  predicted.  This  prediction  is  based  on  the  stress-strain 
curve  for  AGSR  at  70 “F  as  shown  in  Figure  3-12  of  Section  3.2.2.  A  given  element 
in  the  structural  grid  was  assumed  to  have  failed  if  the  level  of  the  hoop  strain 
exceeded  1.35  x  10~^  in/in.  This  is  a  good  approximation  because  the  principal 
strains  in  each  element  are  only  slightly,  larger  than  the  hoop  strains.  Figures 
5-11  and  5-12  indicate  that  the  presence  of  the  groove  at  the  downstream  end  of 
the  Insert  causes  a  localized  region  of  failure. 

5. 6. 1.2. 2  Plastic  Analysis  vs.  Elastic  Analysis 

The  elastic-plastic  analysis  (Case  6)  was  performed  to  determine  how 
plastic  deformation  affected  the  level  of  total  strain  and  displacement  in  the 
Insert  as  a  result  of  the  cooldow~  process.  The  material  yield  properties  and 
tangent  moduli  were  discussed  in  Section  3.2.2,  and  the  method  of  plastic  analy¬ 
sis  was  described  in  Section  5.1.  The  numerical  calculations  indicate  that 
plastic  deformation  reduced  the  stresses  in  the  insert  only  slightly;  however, 
the  strains  and  displacements  were  increased  far  more  significantly.  The  hoop 
strains  as  calculated  by  the  plastic  and  elastic  analyses  are  plotted  as  func¬ 
tions  of  radial  position  in  Figures  5-8,  5-9,  and  5-10.  The  total  strain  in  a 
plastic  analysis  is  defined  as  the  sum  of  the  elastic  and  plastic  portions. 
However,  note  that  the  elastic  portion  of  a  plastic  analysis  is  not  the  saune  as 
the  total  strain  in  an  elastic  analysis.  It  is  seen  from  these  figures  that 
the  maximum  strains  occur  in  the  entrance  plane  (upstream  surface)  of  the  in¬ 
sert,  and  the  lowest  strains  occur  in  the  throat  plane.  The  important  conclu¬ 
sion  to  be  drawn  from  these  figures  is  that  the  total  strains  as  determined  by 


5-9 


the  plastic  analysis  are  significantly  higher  (as  much  as  21  percent)  than 
those  determined  by  an  elastic  analysis. 

Although  the  stresses  determined  by  the  plastic  analysis  are  slightly 
lower  than  those  calculated  by  an  elastic  analysis  (e.g.,  see  Figures  5-5  and 
5-13) f  it  is  the  strain  level  that  causes  material  failure  not  the  stress  level. 
Therefore,  the  fact  that  an  elastic  analysis  is  more  conservative  in  the  estima¬ 
tion  of  stress  levels  does  not  mean  that  an  elastic  analysis  is  more  conserva¬ 
tive  in  the  prediction  of  failure.  Hence  a  plastic  analysis,  which  is  more  con¬ 
servative  in  the  estimation  of  strain  levels  than  an  elastic  analysis,  is  more 
appropriate  for  the  prediction  of  a  material  failure.  The  importance  of  the 
previous  statement  becomes  more  apparent  if  reference  is  made  to  Figure  5-14. 

The  results  of  the  elastic  analysis  indicate  that  failure  occurs  in  a  small 
region  near  the  ends  of  the  insert  and  near  the  AGSR-PG  interface.  This  region 
is  shown  in  Figure  5-12.  The  plastic  analysis  indicates  that  failure  occurs  in 
a  larger  region.  The  failure  criterion  for  the  substrate  was  again  based  on  the 
hoop  strain  state  exceeding  1.35  x  10  ^  in/in. 

Another  difference,  although  relatively  minor,  between  the  plastic  and 
elastic  analyses  is  illustrated  by  Figure  5-15.  This  plot  depicts  the  displace¬ 
ment  of  the  outer  and  inner  surfaces  of  the  AGSK  as  functions  of  axial  position 
for  the  heating  and  cooling  portions  of  the  insert  fabrication  process.  Dis¬ 
placements  due  to  the  cooling  portion  were  determined  by  both  plastic  and  elastic 
analyses.  It  is  seen  that  the  contraction  predicted  by  a  plastic  analysis  is 
slightly  less  than  that  predicted  by  an  elastic  analysis. 

5. 6. 1.3  Machining  the  Ends 

The  final  step  in  the  fabrication  process  consisted  of  machining  0.1  inch 
from  the  upstream  and  downstream  surfaces  of  the  insert  at  70®F.  Since  a  non¬ 
zero  stress  state  existed  in  the  insert  at  70 °F,  the  effect  of  machining  on 
this  state  was  considered.  This  was  done  by  removing  elements  from  the  struc¬ 
tural  grids  (with  and  without  grooves)  to  simulate  the  machining  process  and 
allowing  the  stresses  in  the  insert  to  redistribute.  The  results  from  this 
analysis  indicated  that  the  machining  process  caused  a  negligible  redistribution 
of  stresses,  strains,  and  displacements  in  the  insert. 

5.6.2  Motor  Firing 

The  structural  analysis  of  the  motor  firing  was  performed  by  calculating 
the  stress  state  in  the  insert  at  various  instants  in  time  during  the  motor  burn¬ 
ing.  The  times  chosen  were  1.5,  3.0,  and  60.0  seconds. 
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The  1.5  second  stress  state  was  analyzed  because  the  temperature  gradient 
In  the  insert  Is  most  severe  at  the  beginning  of  the  firing.  Hence,  the  PG 
coating  is  at  a  much  higher  temperature  near  the  surface  than  at  the  AGSR-PG 
Interface.  The  30.0  second  condition  was  important  because  rocket  motor  experi¬ 
mental  tests  were  conducted  for  both  30  and  60  second  durations;  hence,  a  nu¬ 
merical  analysis  of  the  stress  state  at  30  seconds  enabled  investigation  of  an 
intermediate  condition  for  one  firing  as  well  as  a  burnout  condition  for  another 
firing.  The  60.0  second  condition  was,  of  course,  important  because  the  tempera¬ 
ture  of  the  insert  was  highest  near  the  end  of  the  firing.  This  temperature  causes 
the  compressure  force  between  the  insert  and  its  backside  boundaries  to  be  a 
meucimum.  In  addition,  the  differential  expansion  between  the  PG  and  AGSR  would 
be  expected  to  cause  a  high  stress  level  at  the  interface. 

Because  the  material  properties  vary  with  temperature,  the  reference  or 
initial  state  of  the  insert  for  the  stress  analysis  had  to  be  stre8.<}  free  which 
meant  that  the  state  at  3634"F  had  to  be  the  reference  state.  However,  the 
structural  results  from  the  plastic  analysis  showed  that  significant  plastic 
strains  occurred  during  cooldown  of  the  insert.  This  means  that  proper  definition 
of  the  motor  firing  strains  requires  the  input  of  plastic  strains  into  the  stress 
analysis  as  initial  conditions  at  an  insert  temperature  of  70 *F.  However,  the 
handling  of  this  initial  condition  is  beyond  the  capability  of  the  DOASIS  code, 
but  it  can  be  handled  by  the  OASIS  code  (Ref.  26)  ,  although  at  a  significant 
increase  in  cost  which  made  it  impractical  for  this  investigation.  As  a  result 
the  influence  of  the  plastic  strains  on  the  motor  firing  strains  were  neglected 
and  the  state  at  3634‘F  was  used  as  the  stress  free  state. 

As  mentioned  in  Section  4.3,  two  different  nozzle  Insert  temperature 
distributions  were  generated  in  the  heat  transfer  analysis.  One  was  determined 
considering  only  the  nozzle  Insert  and  the  other  was  determined  considering  the 
entire  nozzle  assembly.  For  all  but  one  of  the  stress  calculations  (Case  11) , 
the  former  temperature  field  was  used. 

In  the  following  sections,  the  results  of  the  stress  analyses  are  presented 
separately  for  each  point  in  time.  Tables  5-2  and  5-3  present  a  summary  of  the 
various  conditions  for  which  analyses  were  performed. 

5. 6. 2.1  Stress  State  -  1.5  Seconds 

Stress  analyses  were  performed  at  the  1.5  second  state  for  three  differ¬ 
ent  conditions.  Cases  7  and  8  shared  identical  boundary  conditions;  that  is, 
the  upstream  and  downstream  surfaces  v;ere  pressurized,  and  the  outer  surface 
was  not  allowed  to  displace  radially.  However,  Case  7  analyzes  an  insert  with 
grooves,  and  Case  8  analyzes  an  insert  without  grooves.  Case  9  was  identical 
to  Case  7  except  that  the  reference  stress  free  state  was  chosen  as  70*F.  This 
case,  in  fact,  was  the  only  one  in  which  the  reference  stress  free  state  was 
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not  choa—  m  3C34*f  •  TIm  rasulta  of  Caao  9  were  used  as  a  baeie  for  coaiparlng 
the  paroeent  aoletloa  teohnlqiie  with  previous  superposition  solutimi  techniques. 
The  CMMlts  of  this  ocayarison  are  discussed  in  detail  in  Section  5.6.3.  The 
raaialta  Cor  Caass  7  and  •  are  shown  ^n  Figures  5-16  through  5-21  where  the  shear 
strasaas  at  tha  PG-’hOW  iatarfaca  are  plotted  as  functions  of  axial  position  and 
tha  hoop  straaaae  in  tha  PG  and  AGSR  are  plotted  as  functions  of  radial  position 
at  tha  aatraaca*  throat,  and  exit  planes  of  the  insert.  The  naxisnui  shear  stress 
(Pifora  5*16)  is  ssuillar  in  tha  insert  with  grooves.  However,  it  appears  that 
tha  downstraasi  qroova  is  not  as  effective  as  the  upstream  for  reducing  shear 
strasaas. 

Tha  asrisasi  hoop  stresses  in  the  PG  occur  near  the  nozzle  throat  plane. 
Oosyaring  tha  hoop  strasaas  in  the  grooved  and  nongrnovcd  inserts,  it  is  appar¬ 
ent  that  although  the  distributions  of  stress  are  different,  the  swxJjsun  nagnl- 
tade  of  the  stresses  are  not  very  different.  The  same  conclusion  can  be  drawn 
about  tha  hoop  stress  level  in  the  AGSR  (Figures  5-20  and  5-21). 

5.6.2. 2  Stress  State  -  30.0  Seconds 

MUBwrical  calculations  were  performed  at  the  30  second  state  for  six 
different  sets  of  conditions  (Casas  10  through  15) .  Reference  to  Tables  5-2 
and  5-3  provides  a  susmary  of  these  conditions.  The  temperature  field  deter- 
ainad  by  the  heat  transfer  analysis  of  the  nozzle  insert  was  used  for  all  cases 
except  Case  11.  The  tesyerature  field  for  Case  11  was  the  field  calculated  by 
considering  heat  transfer  to  the  entire  nozzle  assembly.  The  upstream  and  down¬ 
stream  surfaces  of  the  insert  were  always  assumed  pressurized,  except  for  Case 
12  where  a  nixed  boundary  condition  was  placed  on  the  downstream  surface.  From 
this  case,  the  portion  of  the  surface  between  the  outer  edge  and  the  groove  was 
not  allowed  to  displace  axially.  This  condition  was  described  in  detail  in 
Section  5.4.4.  The  outer  surface  of  the  insert  was  not  allowed  to  displace 
radially  for  Cases  10,  11,  and  12.  However,  in  Cases  13  and  14,  the  outer  sur¬ 
face  was  allowed  to  displace  radially  outward  a  distance  of  5  mils  and  16  mils 
^sspectively  before  encountering  a  rigid  boundary,  in  Case  15,  a  0.2  inch  gap 
filled  with  grafoil  was  placed  between  the  outer  surface  of  the  insert  and  the 
silica  phenolic.  The  silica  phenolic  was  assumed  to  be  rigid  with  respect  to 
the  radial  and  axial  displacements. 

The  structural  finite  element  grid  of  the  insert  for  the  30  second  state 
was  slightly  sndified  to  account  for  the  material  recession  at  the  PG  surface. 

A  layer  of  finite  eleswnts  10  mils  thick  was  removed  from  the  structural  grid 
to  approx ism te  the  amterial  recession.  Actual  calculations  of  SMterial  reces¬ 
sion  ware  discussed  in  Section  4.3.2.  Case  10  was  considered  as  the  "baseline” 
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condition  for  the  30  second  state  calculations.  Case  11  was  performed  to  assess 
the  difference  in  the  stress  state  using  the  nozzle  assembly  temperature  distri¬ 
bution.  Cases  12,  13,  14,  and  15  were  performed  to  determine  the  effects  of 
various  boundary  conditions  on  the  stress  state  in  the  insert. 

The  results  of  all  these  calculations  are  shown  in  Figures  5-17  through 
5-23.  They  indicate  that  the  nozzle  assembly  temperature  distribution  produces 
hoop  and  radial  stresses  significantly  higher  than  those  produced  by  the  nozzle 
Insert  temperature  distribution.  The  increa  Ad  stresses  are  principally  due  to 
the  additional  heat  flux  through  the  ends  of  the  insert  from  the  rest  of  the 
nozzle  assembly. 

Comparisons  of  various  outer  surface  boundary  conditions  (Figure  5-22) 
indicate  that  the  highest  compressive  hoop  stresses  (9300  psi)  occur  in  the  AGSR 
if  the  outer  surface  is  not  allowed  to  expand  radially.  The  minimum  stresses 
(-450,  -flOO  psi)  occur  in  the  AGSR  if  the  outer  surface  is  allowed  to  expand 
radially  outward  5  mils.  The  maximum  tensile  hoop  stresses  occur  in  the  AGSR 
if  the  outer  surface  is  allowed  to  expand  radially  outward  16  mils,  or  if  the 
0.2  inch  thick  layer  of  grafoil  is  placed  in  the  gap  between  the  outer  surface 
and  the  silica  phenolic.  The  hoop  stresses  in  the  P3  (Figure  5-23)  were  maxi¬ 
mized  in  compression  if  the  c’.ter  surface  was  not  allowed  to  expand  radially. 

The  hoop  stresses  became  smaller  as  the  outer  surface  was  allowed  to  expand 
further.  Comparison  of  the  results  from  the  two  cases  which  considered  the 
variation  in  the  downstream  boundary  condition  showed  that  this  variation  had 
a  negligible  influence  on  the  stress  distribution. 

In  general,  the  stress  level  at  the  30  second  state  is  higher  than  the 
level  at  the  1.5  second  state.  The  stresses  increase  the  most  at  the  PG-AGSR 
interface.  The  only  stresses  which  decreased  below  their  level  at  the  1.5  sec¬ 
ond  state  were  the  interlaminar  shear  stresses  in  both  the  PG  and  AGSR. 

3. 6. 2. 3  stress  State  -  60.0  Seconds 

Numerical  calculations  were  performed  at  the  60  second  state  for  three 
different  sets  of  conditions  (Cases  16  through  18) .  Reference  to  Tables  5-2  and 
5-3  provides  a  summary  of  these  conditions.  The  only  differences  between  the 
three  conditions  were  the  insert  outer  surface  boundary  conditions.  For  Case  16, 
the  outer  surface  was  not  allowed  to  displace  radially.  In  Cases  17  and  18,  the 
outer  surface  was  allowed  to  displace  outward  5  mils  and  16  mils  respectively. 
Another  10  mils  were  removed  from  the  structural  grid  for  the  60  second  state  to 
account  for  further  material  ablation.  The  total  amount  of  material  removed 
since  the  beginning  of  the  firing  was  20  mils. 

The  stress  state  for  each  of  the  three  60  second  conditions  is  sho%m  in 
Figures  5-17  through  5-21  and  Figures  5-24  and  5-25.  The  latter  two  figures 
provide  a  comparison  of  the  hoop  stress  level  in  the  insert  for  various  outer 
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surface  boundary  conditions.  It  is  apparent  that  the  stress  level  is  the  lowest 
if  the  outer  surface  is  allowed  to  expand  radially  16  mils  before  encountering 
a  rigid  boundary,  and  highest  if  it  is  not  allowed  to  expand  radially  at  all. 

The  hoop  stresses  at  the  PG-AGSR  interface  reach  a  level  at  the  60  sec¬ 
ond  state  which  is  higher  than  at  any  other  time  in  the  motor  firing.  The  peak 
hoop  stress  encountered  for  any  boundary  condition  was  -15200  psi. 

5.6.3  Comparison  of  Present  Technique  with  Superposition  Principle 

Some  previous  stress  analyses  conducted  by  other  investigators  have  been 
performed  on  the  rocket  nozzle  insert  using  a  superposition  principle.  These 
analyses  calculated  the  motor  firing  stresses  by  adding  the  results  of  two 
separate  calculations.  The  first  calculation  generated  the  stresses  resulting 
from  the  cooldown  process  assuming  that  3634 “F  was  a  stress  free  state.  The 
second  calculation  generated  stresses  resulting  from  exposing  an  insert  to  a 
motor  firing  condition  assuming  that  70*^  was  a  stress  free  state. 

The  stress  analysis  technique  used  for  the  calculations  in  this  report 
assumed  that  the  stress  free  state  of  the  insert  was  at  3634 °F  for  both  the  cool¬ 
down  and  motor  firing  condition.  The  validity  of  this  approach  is  based  on  two 
assumptions.  The  first  assumption  is  that  machining  of  the  ends  of  the  insert 
prior  to  placing  it  in  the  nozzle  assembly  does  not  appreciably  alter  the  stress 
distribution  in  the  insert.  The  results  of  Case  4  verify  that  this  assumption 
is  well  founded.  The  second  assumption  is  that  the  effect  of  the  cooldown 
plastic  strains  can  be  neglected. 

If  the  material  properties  of  the  PG  and  AGSR  were  not  temperature  depen¬ 
dent  and  if  the  geometry  of  the  insert  were  not  different  at  the  70*F  and  3634®F 
states,  the  two  techniques  would  produce  identical  results.  However,  this  is 
not  the  case  because  the  partial  differential  equation  governing  the  state  of 
stress  in  the  insert  is  dependent  on  material  properties.  Superposition  of 
solutions  of  this  partial  differential  eqxiation  is  valid  only  if  the  superposed 
solutions  are  governed  by  the  same  partial  differential  equation.  As  shown  in 
Section  3.2.2,  the  material  properties  of  PG  and  AGSR  vary  considerably  with 
temperature.  Hence,  the  partial  differential  equations  governing  each  step  in 
the  superposition  technique  are  different.  Therefore,  it  is  not  valid  to  super¬ 
pose  the  solution  of  each  step.  A  secondary  reason  why  superposition  is  invalid 
is  that  the  insert  contracts  in  going  from  3634 ®F  to  70®F;  hence,  the  region 
governed  by  the  partial  differential  equation  changes  in  cabh  step  of  the  super¬ 
position  technique. 

A  comparison  of  the  results  obtained  by  each  method  is  shown  in  Figure 
5-26.  The  hoop  stresses  in  the  PG  at  the  1.5  second  state  of  the  motor  firing 
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are  plotted  as  functions  of  the  distance  normal  to  the  surface  of  the  PG  using 
both  techniques.  The  results  Indicate  that  the  superposition  principle  diverges 
from  the  present  technique  In  the  critical  region  of  largest  temperature  gradi¬ 
ent.  At  the  PG  surface,  the  results  differ  by  35  percent. 
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TABLE  5-2 


SUMMARY  OF  STRUCTURAL  ANALYEEL 
PERFORMED  USING  DOASIS  COMPUTER  PROGRAM 


Run 

No. 

Solution 

Type 

Substrate 

Time 

Temperature 

Motor 

Firing 

Time 

(bcc) 

Boundary 
Conditions 
(See  Table 
5-3) 

Zero 

Stress 

State 

Remarks 

1 

Elastic 

Grooved 

3634«F 

— 

Stress  Free 
Boundaries 

70*F 

Define  stresses  and  thermal 
expansion  for  substrate  only  ' 

2 

Elastic 

Non- 

Grooved 

3634*r 

Stress  Free 
Boundaries 

70*F 

Same  as  for  Run  No.  1 

1 

3 

Elastic 

Grooved 

70*F 

Stress  Free 
Boundaries 

3634*F 

Cooldown  following  application  i 
of  PG  coating 

4 

Clastic 

Grooved 

70*F 

Stress  Free 
Boundarie.>: 

3634*r 

Define  stress  redistribution 
due  to  machining  ends 

5 

Elastic 

Non- 

Grooved 

70*F 

— 

Stress  Free 
Boundaries 

3634»F 

Same  as  for  Run  No.  3  ^ 

6 

Plastic 

Grooved 

70*F 

Stress  Free 
Boundaries 

3634*F 

Plastic  cooldown  analysis 

7 

Elastic 

Grooved 

Nozzle 

Insert 

1.5 

Tl,  0), 

I'l,  D) 

3634*F 

Motor  firing  stress  analysis 
with  zero  backside  disp. 

6 

Elastic 

Non- 

Grooved 

Nozzle 

Insert 

1.5 

n,  01, 

Ifl  ,  01 

3634*F 

Same  as  for  Run  No.  7  j 

9 

Elastic 

Grooved 

Nozzle 

Insert 

1.5 

n,  01, 

Ul  ,  01 

70»F 

Same  as  for  Run  No.  7  I 

70*P  assumed  as  stress  free  state] 

10 

11 

Clastic 

Elastic 

Grooved 

Grooved 

Nozzle 

Insert 

Nozzle 

Assembly 

30 

30 

12,  01, 
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SECTION  6 


EVALUATION  OP  THERMOSTRUCTURAL  RESULTS 


The  Structural  results  presented  in  Section  5  show  that  the  primary  struc¬ 
tural  problems  are; 

•  The  interleuninar  shear  stress  at  the  interface  between  the  AGSR  sub¬ 
strate  and  the  PG  coating  during  cooldown  of  the  insert  during 
fabrication 

•  The  tensile  hoop  stress  or  strain  in  the  AGSR  during  cooldown  of  the 
substrate  during  fabrication 

•  The  compressive  hoop  stress  or  strain  in  the  PG  coating  during  the 
motor  firing 

•  The  definition  of  the  constraint  on  the  outside  diameter  of  the  AGSR 
substrate  in  order  to  minimize  the  hoop  stress  in  the  PG  coating 
while  at  the  same  time  keeping  the  stresses  in  the  substrate  to  an 
acceptable  level. 

The  primary  problem  in  evaluating  the  structural  results  is  that  the  structural 
analytical  model  used  in  the  DOASIS  computer  code  has  not  been  validated  for  the 
PC  coated  nozzle  insert  due  to  the  scarcity  of  meaningful  experimental  data  for 
comparison.  The  "pass  or  fail"  type  of  data  which  results  from  the  motor  firing 
indicates  only  whether  or  not  the  structural  limitations  of  the  insert  have  been 
exceeded.  Since  these  data  very  seldom  provide  any  information  concerning  the 
margin  of  safety  of  the  insert  if  the  test  was  successful  or  how  severely  the 
insert  was  overstressed  if  the  test  was  a  failure,  the  data  only  provide  an  end 
check  on  the  validity  of  the  structural  analysis.  In  addition,  since  the  struc¬ 
tural  analysis  is  based  both  on  the  material  mechanical  and  thermal  properties 
and  on  the  structural  analytical  model,  the  source  of  analysis  error  can  not  be 
identified  solely  on  the  basis  of  "pass  or  fail"  motor  firing  data.  To  identify 
the  source  of  error,  quantitative  data  at  various  stages  during  the  fabrication 
process  are  required  for  comparison  with  analytical  predictions.  The  only  data 
of  this  type  which  is  currently  available  for  the  insert  analyzed  in  Figure  5-2 
are  outside  diameter  measurements  of  the  AGSR  substrate  prior  to  and  after  the 
deposition  of  PG  coating.  The  experimental  net  displacements  were  40,  20,  and 
17  mils  for  the  upstream,  throat,  and  downstream  insert  locations.  The  trend  of 
these  measurements  are  in  agreement  with  the  predicted  displacements  presented 
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in  Figure  5-15,  but  the  predicted  displacements  are  a  factor  of  approximately 
four  low.  However,  it  should  be  noted  that  the  high  value  of  net  displacement 
occurs  at  the  upstream  end  of  the  insert  which  is  also  the  region  where  extensive 
failure  was  predicted  by  the  plastic  analysis  (refer  Figure  5-14) .  A  fact  of 
equal  significance  is  that  if  the  measured  displacements  mentioned  above  are 
correct,  then  the  maximum  tensile  hoop  strain  at  the  outside  diameter  of  the 
substrate  is  approximately  11.5  mils.  This  value  of  strain  is  a  factor  of  approx¬ 
imately  10  greater  than  the  strain  at  fracture  based  on  uniaxial  measurements 
(refer  Figure  3-12) . 

The  above  results  become  more  difficult  to  interpret  based  on  recent 
thermal  expansion  data  which  were  obtained  from  Reference  27  and  presented  also 
in  Figures  6-1  and  6-2.  The  data  presented  in  Reference  27  was  obtained  by 
Southern  Research  Institute  on  PG  samples  fabricated  by  Atlantic  Research  Corpor¬ 
ation  and  thus  should  be  representative  of  the  PG  in  the  nozzle  throat  insert. 

The  AGSR  data  are  also  based  on  material  samples  typical  of  that  used  in  the 
nozzle  Insert.  The  data  shown  in  Figure  6-1  compare  the  recent  "a"  direction 
PG  and  AGSR  thermal  expansion  data  with  that  presented  previously  in  Figure  3-19. 
The  recent  "c"  direction  PG  thermal  expansion  data  are  compared  in  Figure  6-2  with 
that  presented  previously  in  Figure  3-20.  The  data  presented  in  Figures  6-1  and 
6-2  show  that  the  recent  "c"  direction  PG  thermal  expansion  data  are  in  good 
agreement  with  previous  data  but  that  the  recent  "a"  direction  PG  thermal  expan¬ 
sion  data  are  in  very  poor  agreement  with  previous  data.  If  indeed  the  data  on 
ARC  PG  presented  in  Figure  6-1  are  representative  of  the  thermal  expansion  in 
the  PG  coated  nozzle  insert  then 

•  The  mismatch  of^the  PG  coating  and  the  AGSR  substrate  thermal  expan¬ 
sions  is  minimal  below  the  deposition  temperature  of  3654®F  and,  thus 
the  cooldown  stresses  should  be  minimal 

•  The  compressive  stresses  in  the  PG  coating  during  the  motor  firing 
would  be  much  larger  than  those  presented  in  Section  5. 

The  discussion  presented  in  the  previous  paragraphs  does  not  shed  much 
light  on  interpreting  the  results  presented  in  Section  5.  However,  if  the 
assumption  was  made  that  the  property  data  used  in  the  analysis  were  correct, 
then  the  following  conclusions  can  be  made 

•  The  predicted  cooldown  hoop  strains  in  the  AGSR  indicate  structural 
failure  of  the  substrate. 

9  The  PG  coating  would  be  structurally  adequate  if  a  16  mil  expansion  of 
the  substrate  was  allowed  during  the  motor  firing.  However,  the  pre¬ 
dicted  hoop  stresses  in  the  PG  coating  indicate  structural  failure  if 
the  allowable  expansion  is  on  the  order  of  5  mils. 
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